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Introduction 



Inflation [1, 2] stands as the central paradigm of modern cosmology. It earned its 
place by automatically solving the so-called flatness, horizon and monopoles problems 
of standard Big-Bang cosmology. It further provides a mechanism for the generation 
of density perturbations in the early Universe, which lead to Large Scale Structures 
(LSS) [3]- [7] in the distribution of galaxies and temperature anisotropies in the Cosmic 
Microwave Background (CMB) [8]- [12]. 

Over the years, many inflationary models have been put forward that account 
for the spectrum of primordial perturbations. These perturbations are to a good 
approximation Gaussian as they are well-described by the leading quadratic terms 
in the inflationary Lagrangian. On the other hand, studying higher order interaction 
terms is crucial in characterizing the various inflationary mechanisms and in pursuing a 
deeper understanding of the inflationary dynamics. With the advent of new generation 
experiments (the launch of the Planck sateUite [13, 14] and the continued analysis of 
WMAP data [15]), which provide improved sensitivity to deviations from Gaussian 
statistics, such a theoretical investigation is being timely complemented by observations 
and is rendered all the more urgent and important. 

In studying primordial non-Gaussianitcs [16, 17, 18, 19] one aims at handling quantities 
that arc conveniently related to obscrvablcs and that, at the same time, have an 
immediate expression in terms of higher order operators in the inflationary action. This 
is the case for n-point functions of the (scalar) degree of freedom that drives inflation. 
Starting with the three-point function, these correlators give an explicit measure of 
the deviation from Gaussianity and depend directly on the interactions in that they 
would vanish if we were to truncate the Lagrangian at quadratic order in perturbations. 
Althought in this work wc will be concerned only with higher order correlators at tree 
level, it is important to point out that interactions give also a (small) correction to the 
power spectrum of curvature perturbations once loop calculations [20] are considered, 
so that too is a measure of non- Gaussianity. 

The process of characterizing the various inflationary mechanisms according to their 
non-Gaussian properties can be schematically outlined. First, one determines if the 
amplitude of the n-point function, or, more conveniently, its Fourier transform, is large 
enough to be detectable. In the affirmative case, it is instructive to proceed and study 
the complete dependence of the correlator on the external momenta, the so called 
shape-function. It is in fact crucial to know in what type of momenta configuration 
the shape-function peaks. Indeed, in analyzing the data one must know beforehand 
if a specific shape-function is expected so that the appropriate corresponding non- 
Gaussianity estimator can be built. Not doing so would result in a system which is 
essentially blind to a whole family of inflationary models (or, in an effective fleld theory 
approach, a family of interaction operators) which might well be there and be important. 
In this work we aim at differentiating among the many inflationary mechanisms trough 
the study of their non-Gaussianities. As a consequence, it seems natural to look for an 
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approach that captures all the general features of an inflationary theory and reduces to 
specific models in appropriate regimes. We specify here that we are limiting our analysis 
to the case of inflation driven by a single scalar degree of freedom and therefore will 
not discuss vectors as has recently been done, for example, in the context of statistical 
anisotropies. As one might suspect, an approach that allows for a unifying perspective 
on inflation can be found in the realms of effective fleld theories. 

A precise prescription based on these ideas has been put forward in [21], whose approach 
we will closely follow in this work. In [21] the authors give, subjected to mild caveats, 
the most general action for inflation driven by a single scalar degree of freedom. The 
Lagrangian for this theory is written down in detail and it turns out one is able to re- 
obtain most specific infiationary models by switching on or off appropriate coefficients 
driving various interaction operators in the action. The hope is to be able to pin 
down the specific interaction terms generating distinctive features in the bispcctrum and 
trispectrum of curvature perturbations. This could result in observational bounds on the 
value of the coefficients (we generically call them Mn) driving the various interactions 
at third and fourth order I. The power of the effective field theory approach is quite 
manifest in that, in principle, it allows these coefficients considerable more freedom 
than what they are granted in any specific inflationary model. In fact, by being for the 
most part free parameters (a couple of these coefficients arc to obey some inequalities 
if one wants, as we do, the generalized speed of sound to be smaller than unity §), 
the Mn' s allow for the description of known interactions with relative weights which 
would otherwise be fixed. By employing effective field theory then one enlarges the 
region of the parameters space than can be spanned. Besides that, in the complete 
effective Lagrangian some of the Mn coefficients multiply {extrinsic curvature- generated) 
operators that are sometimes neglected in the literature but should in principle be 
studied as, in fact, their contribution to higher order correlators can be relevant and 
this effectively increases the dimension of the parameters space of the theory. 
Quite interestingly, the formalism of [21] also sheds new light on effects due to 
symmetries in the action: for example, one can clearly see how a reduced speed of sound 
often automatically results in an enhanced non-Gaussianity. One more advantage that 
comes with employing the proposed setup is of calculational nature: in the so called 
decoupling regime (which implies working in a specific energy range) the dynamics of 
the metric decouples from the one of the scalar that drives inflation thus rendering the 
Lagrangian itself and the higher order correlators much easier to handle and calculate. 
This mechanism is very reminiscent of what happens in standard quantum fleld theory 
and goes under the name of equivalence theorem. 

In this work we aim to integrate and summarize the investigations presented in 

I On the other hand, as we will see, this approach proves itself useful already at second order in 

pcrtiirbations. 

§ One should also mention the general requirement that stems from working within the regime of 
validity of perturbation theory: namely the Mn coefficients are to be smaller than the mass of the 
underlying theory. 
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[22, 23, 24]. The common ground these papers share rehes on the use of effective field 
theory methods within a general effort to characterize as many inflationary models as 
possible according to the non-Gaussian properties they exhibit. These properties take 
the form of various predictions for cosmological observables such as the power spectrum, 
its tilt, its running; the bispectrum amplitude, Jnl, its running; the trispectrum 
amplitude. 

As mentioned, already at second order in perturbations the effective theory approach 
reveals interesting facts. In [22, 24], making full use of the freedom on the Mn's, we are 
able to write a very general quadratic Lagrangian which captures the quadratic theory of 
the entire class of the so called generalized slow-roll theories (often referred to as P{X, (j))- 
model) and also covers models of inflation for which terms with more than one derivative 
acting on the scalar driving inflation are allowed (e.g. Ghost inflation [25]). From the 
general quadratic action we move to the equation of motion: we solve it by imposing 
the Bunch-Davies vacuum condition and by requiring to re-obtain the known solutions 
for specific infiationary mechanisms in the corresponding limits. The resulting power 
spectrum is now a function of five, rather than just three, independent coefficients. This 
in turn means that its tilt depends on the usual three slow- roll parameters but also on two 
additional coefficients, two generafized slow-roll parameters. These results are presented 
at leading order and, in particular cases, at next-to-leading order in generalized slow-roll 
approximation. 

We then turn our attention to non-Gaussianities (NG), starting with the bispectrum. It 
is important to stress here that we will consistently study mostly models that generate 
large, detectable NG and, in doing so, it will become clear that there is an important 
region of the parameters space of the effective theory where the contribution of (often 
neglected) extrinsic curvature-generated || interaction terms is important and, possibly, 
leading. We start from the most general third-order effective action that originates 
from the prescription of [21]. Employing the IN-IN formalism [26, 28, 27, 29], a 
thorough analysis of the various terms contribution to the amplitude, /atl, and shape of 
the three-point correlator is performed. One immediately reproduces results of the 
current literature and re-discovers known features such as the fact that large non- 
Gaussianities may be generated if the sound speed (now a generalization thereof) is 
much smaller than unity. On the other hand, a number of noteworthy novelties arise 
which are due precisely to the effect of extrinsic curvature-generated terms. Upon 
requiring a small generalized speed of sound and using at full the freedom on those 
Mil coefficients that first appear at third order in perturbations, one finds that the 
leading contribution to f^L can indeed come from curvature contributions and there's 
more: the shape-function that some of these terms generate has peculiar, distinctive 
features. In P{X, 0)-models of inflation, but also in Ghost inflation, the typical shape- 
function generated by a single leading interaction term in the cubic Lagrangian will peak 
in the so called equilateral conflguration (that is, in Fourier space, all three external 

II In the approach of [21] the gauge choice identifies a preferred slicing of spacetime; the extrinsic 
curvature tensor K^^, describes the geometry of such slicings. 
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momenta are equal ki = k2 = k^). Interestingly, a number of curvature terms we 
study generate a shape-function that peaks in the flat configuration {ki = 2/c2 = 2ks). 
Such a configuration is quite uncommon for single-field models of inflation and cannot 
be obtained without considering extrinsic curvature-generated terms in the action 
unless one relaxes the Bunch-Davies vacuum condition on the wavefunction [30] or 
considers linear combinations of third-order operators [31]. Our study then enlarges 
the classes of single-field inflationary models whose bispectrum signature may consist 
in a flat shape-function and makes this feature somewhat more natural as a flat shape 
now independently originates from several interaction terms within the B-D vacuum 
condition for the wavefunction. 

Prompted by the /jVL-related findings concerning extrinsic-curvature interactions, we 
investigated the corresponding contribution to the running of J'nl itself. This quantity 
is expected to be small, roughly of the order of the generalized slow-roll parameters. 
The actual calculation shows that the contribution from curvature terms might be the 
leading one for the running and that it can be larger than the slow-roll parameters since 
here the quantities involved are both (quite constrained) Mn coefficients which appear 
in the quadratic Lagrangian and Mn's that first show up in the cubic action. 
Following [23], we present what follows quite naturally after the investigations briefly 
outlined above, namely a study of the trispectrum generated by the effective action up 
to fourth order. Here again, the novelties are to be found in the effects of the extrinsic 
curvature-generated interaction terms. These can provide the dominant contribution 
to the four-point function and present very distinctive patterns in the form of the 
trispectrum shape-function. 

The powerful and very convenient setup that we are going to thoroughly describe below 
is subjected to some specific hmitations. We will discuss them as they arise along the 
presentation. 

The paper is organized as follows. In the first section we briefiy review the setup 
of [21]. In Section 2 we proceed with the analysis of the quadratic Lagrangian for 
the perturbations. We solve the corresponding equation of motion and calculate the 
resulting spectrum, its tilt and running. In Section 3 we take on cubic interactions 
obtaining the amplitudes and shape-functions for the leading contributions to the 
bispectrum of curvature perturbations. We also calculate the running of /at/, in a number 
of particularly interesting cases. Section 4 is dedicated to the trispectrum: here most 
of the possible quartic interactions are considered and special attention is devoted to 
operators that are invariant under two specific symmetries recently introduced in the 
literature. In the Conclusions we summarize the significance of our results and comment 
on further work. In the various Appendices we present some explicit calculations which 
have been omitted from the main text for the sake of simplicity. 
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2. The effective action for single-field infiation (up to third order) 



Our goal will eventually be to write down the complete theory of single-field models 
of inflation up to fourth order in perturbations. We will follow the effective theory 
approach first introduced in Ref. [21] of which we now give an outline. 
The single scalar field (f) which one assumes to be responsible for infiation is splitted as 
usual in an unperturbed part, the background, plus a fiuctuating one: 

(l>(x,t)^Mt) + H(x,t). (1) 

At this stage it is essential to underline the gauge choice that is made as it is not the 
most common one found in the literature. In [21] and in the following one works in 
the comoving (or unitary) gauge for which Scj) — (see also [34]), the scalar degree of 
freedom is now hidden in the metric. Once this choice is implemented, the Lagrangian 
will no more be invariant under full spacetime diffeomorphisms (diffs) but only under 
spatial reparamctrizations. This is the starting point to write the most general space 
diffs invariant Lagrangian at the dcsidcrcd order in perturbation theory in an effective 
theory approach. In [21] the authors prove that, once an approximate shift-simmetry is 
required, their second and third order action is the most general one (see also [49] for 
an interesting perspective on the most general effective Lagrangian for inflation). One 
can then use the so called Stueckelberg trick to restore full spacetime reparametrization 
invariance. As a by-product of this latter procedure, the degree of freedom hidden in 
the metric shows up again as a scalar field. 

Let us start from the general theory before full spacetime reparametrization 
invariance has been restored: the most general space diffs-invariant action in unitary 
gauge can be schematically written as [21]: 

J d^x^g F{R^,p,, K^,, V^, t), (2) 

where K^i, is the extrinsic curvature tensor on which we will soon elaborate more and 
the "0" components of the metric tensor g^'^ are free indices. Considering fluctuations 
around a FRW background amounts to studying the following action: 

J d^x^f^ \mI^R + Ml^Hg^ - M|i (^3//^ + if) + 

F^'^Kg'' + 1, 5R,,p,; V,- 1)] , (3) 

n>2 

where the F*^") functions contain fluctuations which are at least quadratic. The next 
step is to restore full spacetime reparametrization invariance. To see how it works, we 
borrow a simple example from [21] and consider the following sample action terms 

j d'x^g[A{t) + B{t)g'\x)\ . (4) 

We are interested in time reparametrization t ^ t + ^^{x^t)] x ^ x, under which the 
above action (after a simple variable redefinition) reads 



j d^x \/-g{x) 



A{t - i^{x)) + B{t - i^{x)) — — g'' {x) 
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(5) 



Upon promoting to a field, ^0(2^) = — 7r(x) and requiring the following gauge 
transformation rule 7r{x) — )■ 7r{x) — ^o{x) on tt, the above action is invariant under 
full spacetime diffeomorphisms. The scalar degree of freedom tt makes its appearance 
in the time dependence of the A, B coefficients and in the transformed metric. More 
into details, under time reparametrization the metric g^'^ transforms as follows: 

which implies 

g^'^g^'; g'' ^ {I + ix)g'' + g'^d.n; (7) 
/° ^ (1 + 7r)V° + 2(1 + 7r)^°'9i7r + g'^dindjU. 

This procedure has been borrowed, conceptually unchanged, from standard gauge 
theory: a Goldstonc boson which transforms non linearly under the gauge 
transformation provides the longitudinal component of a massive gauge boson. At 
sufficiently high energy such Goldstone boson becomes the only relevant degree of 
freedom. This is the so-called equivalence theorem. The same is true for our case: 
for sufficiently high energy the mixing with gravity becomes irrelevant and the scalar 
TT becomes the only relevant mode in the dynamics. This is the so-called decoupling 
regime. Let us clarify this concept with a simple example. Consider the following 
contribution, taken from Eq. (3) 

M|iij^o° ^ M|iij((l + 7r)2^°° + 2(1 + 7r)y°^9,7r + g'^ d.T^dji,). (8) 

We focus on the quadratic part of the first term in the above equation. Upon canonical 
normalization, tTc = Mpiij^/^vr and g^^ — Mpig'^^, one gets 

M|ii7(^o° + 27r^°o + TT^g^') = tt^ + 2H'/''7r,g^' + Mp.Hgf. (9) 

Consider the second term of Eq. (9) which mixes gravity with the scalar. Since tTc ~ -EtTc, 
at energies higher than ~ H^/'^ the term tt^ dominates the dynamics. This turns out to 
be true in general: the number of derivatives (which in Fourier mode would basically 
give an energy-dependent coefficient in front of tt) is higher in terms containing only tt's 
than in the mixed terms and therefore there exists an energy threshold above which the 
scalar decouples from gravity. Since in explicitating the F'^'^^ term in Eq. (3) there can 
be, in principle, other quadratic terms that go like tt^, one has to consider which one 
is the leading kinetic term and determine the canonically normalized field tTc and the 
energy threshold accordingly. To take the safe route, one might well take the energy 
threshold, i?mixj to be the highest one of this set. Since one is concerned with correlators 
just after horizon crossing, one concludes that the decoupling procedure works as long 
as the decoupling energy is smaller than the Hubble rate H . More precisely, we can 
anticipate that the kinetic terms in F^"^^ which are going to matter in our discussion 
come with coefficients M| and Mp^eH^. The condition E'mix < H is then satisfied if 
Mp^eH^ > M|, where e = —H/H'^ is a slow-roll parameter; if this is not the case we 
need to assume M| < Mp^H'^. 
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Prom now on we will work in the decoupling regime. In considering the terms of Eq. (3), 
we will therefore use only the unperturbed entries of the metric tensor. In order to write 
the effective Lagrangian up to, say, third order ^, we start from Eq. (3) and follow the 
algorithm given in [21]. Fluctuations are encoded in the F*^") terms. In order to be as 
general as possible, we also include all possible contributions up to third order coming 
from extrinsic curvature K^y terms. In fact, it is instructive at this stage to step back 
and consider the action in Eq. (2). Given a theory which is space diffs- invariant, one 
can always identify a slicing of spacetime, described by a timelike function t{x), which 
realizes time diffeomorphism: on surfaces of constant t the time symmetry breaking 
scalar is also constant. Before selecting a gauge, there is still the freedom to make a 
choice on t{x) and working in the unitary gauge amounts to requiring t — t. In order 
to describe the geometry of this preferred slicing, one employs the extrinsic curvature 
tensor. In writing down such a tensor, one needs two ingredients: the unit normal 
vector n^, perpendicular to the constant t surfaces, and the induced metric hn^. These 
are defined as 

n^, = —=2^==^ -T==; hi^„ = Qi^y + n^riy, (10) 
which allows us to write 

/Y^^_Ai^v^n^ - 2^_^oo)3/2 + 2(-^oo)5/2 2(-^oo)i/2 ■ 

The above expressions can be used to write explicitly the most generic third order action 
for the fluctuations around the FRW background: 



/ 



+ l^M.itng^' + 1)3 - ^^(/° + 1)SK^, - ^^SK%' 



SRi^JK^SKP^ . (12) 



3! ' ^ 3! 3! 

_ Mg{t ) 
3! 

The coefficients M2,Ms and Mi, • • • ,Mg, to which we will often refer to as the Mn 
coefficients, are in principle generic; we will comment on their physical significance as 
we discuss them more in detail. All the M coefficients multiply extrinsic curvature- 
generated interactions. A given particular set of values (or bounds) for the Mn's will 
specify a given inflationary theory. 

The action, as written in (12), is not yet invariant under full diffeormophisms. One 

^ We will write the explicit expression up to fourth order when concerned with the trispectrum, here 
we limit ourselves to third order for simplicity. 
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needs to follow exactly the steps illustrated in Eqs (4), (5) and (8) and promote to 
a field tt with the proper gauge transformation. 
In the decoupling limit we find: 



\ ) + ^-^(^--^^ 

mitf f idf7r){d]7r) + 2H{d,'KfdjTi + 2'kdl'Kd,Ti\ M-j{t) {dfirf 



2 V «^ / 3! a6 



(13) 



3! ^ ' 3! a6 
A few clarifying comments are in order: 

• If we consider terms only up to second-order, for M2 = Mi,2,3 = one recovers the 
usual quadratic Lagrangian for the fluctuations, with sound speed = 1 and the 
standard solution to the equations of motion. Switching on M2 amounts to allowing 
models with sound speed smaller than unity, 1/c^ = 1 — 2M|/(M|iij), which are 
often hnked to a high level of primordial non-Gaussianity [51, 21]. Furthermore, 
turning on M2,3 in the de Sitter limit, one recovers Ghost inflation [25]. On the same 
lines, keeping all the M's vanishing, but going to third and higher order with the 
M's, one can retrieve the interactions that describe DBI inflation [50, 51, 53]. The 
list of correspondences continues with K-infiation theories and others, thus showing 
how the effective action approach provides a unifying perspective on inflationary 
models [21]. 

• The action in Eq. (13) has already been written with large non-Gaussianities in 
mind. This means that, at every order in fluctuations and for each M and M 
coefficients, we have selected those leading terms which will eventually generate 
large three-point correlators. To clarify this point, we provide a simple example. 
Let us consider the terms up to second order in Eq. (13) and set conveniently 
-^1,2,3 = 0. The properly normalized solution to the equation of motion will be 
the usual 7rfc(r) oc iHe~^^^^'^{l + ikcgr). It is straightforward to verify that, at 
the horizon crossing , tt ~ Hn and Vtt ~ H/cg tt. Therefore, among the tt terms 
with the same number of derivatives, the ones with the highest number of space 
derivatives dominate in the <C 1 limit. Generalizing these estimates for the 
classical solution (which we will describe below) obtained from the equation of 
motion of our complete action, one selects the terms in Eq. (13). 

• there is also the comparison between same perturbative order but different M terms 
to be made. In the literature, all non zero coefficients in front of the various 
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operators are generically assumed to be of the same order (see for example the 
discussion concerning the orthogonal configuration in Rcf. [31] for an interesting 
perspective). We shall not restrict ourselves to this situation. Note that, were the 
coefficients to be all of the same order, one could already identify the dominant 
operators. For example, consider the third order contributions Mf (9i7r)^9|7r/a*^ ~ 
M^{H^ / c^)n'^] for Cs -C 1 this will be a leading contribution with respect to, 
say, M2{tfH{dlT^){djT^f - M"^ H {H^ cI)tt\ This is due to the fact that in the 
effective Lagrangian every additional derivative comes with a H/M <^ 1 factor 
attached: one is basically doing an H/M expansion where M is roughly the energy 
range of the underlying theory. In the last example we have intentionally picked 
terms with the same power of Cs at the denominator. Let us now look at the 
M7{t){df7r)^ ~ M(iJVc^)7r3 term though ; comparing this contribution with the 
Ml term amounts to comparing with MH'^/cl. We see that for a very small 
speed of sound the M7 contribution may still be relevant. These examples justify 
our strategy of including all the terms in Eq. (13) compatible with Cs <^ 1 (we will 
make more comments on this point in the next Section). 
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3. Classical solution and power spectrum 



We now proceed to solve the equation of motion for the second-order effective Lagrangian 
at leading order in slow-roll: 

A = a' {M^^Hid.^r + 2M,V + M!H^-^ - ^^M^f - ^^.{d^.^f) ■ (14) 
Let us not a few facts about the above expression. 

• Eq. (14) is the most general second-order Lagrangian in unitary gauge provided the 
approximate symmetry of the underlying theory is such that only derivative terms 
of TT appear in the action. 

• In order to make contact with more familiar notation, we stress that to first order 
(which is all we need here) the scalar tt is lineraly related to the dimensionless gauge 
invariant quantity C, the curvature, hy ( — —Htt. 

• In full generality the M coefficients above should be time dependent (we will deal 
with such a case in a specific paragraph). However, if one is only interested 
in performing leading-order calculations, then, due to a generalized slow-roll 
approximation, one can safely consider them as constant. 

• The action in Eq. (14) and, as we have seen, its higher order counterparts, are 
generally written with large non-Gaussianities in mind. This results in a number 
of quadratic operators being left out from the formula above. Let us stress 
already at this stage though that, should one decide to include all these subleading 
contributions in the action, the functional expression of the solution will not change, 
one merely redefines a couple of approximately time-independent coefficients. This 
is due to the fact that the types of operator tt^, (9j7r)^, {dji^Y are already saturated 
at the level of Eq. (14). 

We are now ready to tackle the equation of motion. After the usual change of variable, 

— * — * 

Tr{k,t{T)) — a{T)u{k,T), the equation of motion can be written as: 
2 

u" -u + aoPu + (3ok^T^u = 0, (15) 

where f^o are time independent (again, at leading order) dimensionless coefficients. 
This equation has been written in the context of tilted Ghost Inflation [35] and to our 
knowledge, it has not been solved analitically before Ref. [22], where the analytical 
solution has been briefly introduced and used for the computation of the three-point 
function. Here we discuss in much more details the properties of this solution. At this 
stage one can immediately recognize ckq as the more common and f3o as the constant 
a'^H'^/M'^ first introduced in [25]. The complete expression for the coefficients is: 

° -Mj,iH + 2Mf 2(-M|.;ij + 2M|)' ^ ^ 

so that one reobtains the actual for Mi = 0. Note that one can simply look up the 
e.o.m. solution for DBI-like inflation if = = Mi and Ghost Inflation in the de Sitter 
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limit provided ao = 0. Let us pause here to comment on the possibihty of a negative 
CKo (see also [31]). Such a scenario would result in a region in the /c-space, whenever 

laolA;' > ^oA:V - ^, (17) 

for which the solution to the equation of motion will behave exponentially. Such a 
possibility raises a number of issues we address below. First of all, in order to keep 
control of the negative ckq region of the parameters space of the theory in the ultraviolet, 
one requires that the (positive) Pok'^T^ prevails over the ao contribution before k reaches 
the cutoff scale A. Considering that on the IR side, as we will show, the modes will 
eventually freeze outside the horizon, the case of a negative ao should not in principle be 
disregarded. On the other hand, a lot of care should be exerted because an exponential 
phase of the modes for a sufficiently wide k region could generate values for higher order 
correlators that directly contradict available observational data. 

We could now proceed to solve the complete equation of motion but, equipped with 
just equation (15), we can already make some educated guesses on the behaviour of the 
wavefunction. First of all, the typical oscillatory behaviour deep inside the horizon is 

to be expected in this more general case as well: both aok"^ and fiok^r^ cause wave- 
like behaviour (see Fig 1 below) of the wavefunction, while the (— 2/r^) contribution 
is negligible. This is important in that it tells us the main contribution to correlation 
functions will be coming, as usual, from the horizon-crossing region. Note here that, as 
far as (3q ^ 0, the 'Ghost Inflation' term will eventually lead the oscillation if one goes 
deep enough inside the horizon. 

On the other hand, in the r — > hmit, (— 2/r^) will be leading the dynamics and we 
expect to recover the usual, frozen modes. As is familiar from the DBI-like cases, it is 
convenient to introduce the notion of an effective horizon, placing it where the oscillatory 
behaviour stops being dominant. In formulas: 

aoA;' + /3oA;V,' = 4 ^ ^* = , ^ - (18) 

^* k^ ao + + 8^0 

For ^0 = 0, CKo ~ 1 one recovers k'^r^ ~ 1 at the horizon. 

At this stage we can perform a consistency check and show how one can generalize the 
argument, initially borrowed from DBI-like inflationary models, that in comparing terms 
at the same order in perturbations and with the same overall number of derivatives, 
the ones with the most space derivatives are dominating in the ^ 1 limit. The 
generalization of this argument consists in restricting the parameters space to the 
tto ^ 1 md (3o ^ 1 region. Consider Eq. (15) in Fourier space; in full generality 
one expects Vtt ~ /ctt and tt ~ Hn so what needs to be done is relate k with H at the 
horizon. Using equation (18) and r ~ —l/{aH) one obtains 

* . . (19) 

Since the main contributions to correlators comes from the horizon-crossing region, 
this shows that, for {ao, /3o) <^ 1 we can still identify leading terms in the Lagrangian 
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according to the standard procedure. We will strictly follow this procedure when working 
with the cubic and quartic action for the field tt in all the following sections. On the 
other hand, when dealing with the quadratic action, the calculations are simple enough 
so that we can account for all the terms, not just the leading ones. 



3. 1 . Wavefunction 

Let us verify all this quantitatively. The solution to Eq. (15), being of second order, will 
come with two A;-dcpcndcnt integration constants. We have determined their values by 
requiring to re-obtain the known DBI and Ghost solutions in the corresponding limits. 
The general wavefunction reads: 



Uk{r) = 



1 iao 



-4 -VI' -2'-^^^^^^, 



2V4r 



1 iao 



C2{k) , 



(20) 



Where Q stands for the confluent hypergeometric function and C is the generalized 
Laguerre polynomial. We verified that, properly adjusting the integration constants 
according to 



lap 



^-3/2 



■s/Mj^iS H + 2M| 2V4 r 



3/2 — vft) 

' 4(iao+v^) 



C2{k) = 0, 



(21) 



one obtains, in the appropriate limits [56], the wavefunctions of standard inflation and 
Ghost Inflation [25]. We can now write our solution: 



T^kir) = 



i^Mlem + 2M| V27o^'r(| + 



4ao-4i-v/iSo 



(22) 



where 7 = ckq + and T{x) is the Euler gamma function. 

We note in particular that for — 0, Eq. (22) immediately reduces analitically to 
the Ghost Inflation wavefunction 7rfc(r) = {H^-rfl'^ / ^/2Ml)^H\|^{\^/^k'^T'^) with 
T-L^l^ being the Hankel function of the first kind. On the other hand, one can easily see 
numerically that the DBI solution is recovered in the /3o — > limit. 
To give some intuition on the behaviour of the general, interpolating wavefunction, we 
plot it in several (q;o,/3o) configurations. For overall consistency in the comparisons, in 
all the following pictures we have chosen points in the (ao; /9o)-plane so that the horizon 
crossing always lies at the same point, numerically r* = — v^, and we have plotted the 
wavefunction from well inside the horizon (t = — lOr*) uptOT = 0. 
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Figure 1. On the left, the general wavefunction in the DBI-like, /3o — >■ 0, limit; on the 
right the DBI-like solution itself. We plot the real part and find perfect agreement, 
same holds for the imaginary part. To produce the plot the parameters have been set 
to: ao = 0.1, f3o —i' 0, k = 1, H = 1 and the Planck mass-dependent normalization has 
been neglected. The corresponding plot can be omitted for the Ghost limit since in 
that case we recover the Ghost solution analytically. 




Figure 2. On the left, the Ghost Inflation wavefunction (ao = 0). On the right the 
general interpolating solution calculated ioi oq = 1/2 , (3q — 1/4. 



From these plots we learn several things. First, as argued before, there is a common 
oscillatory behaviour once inside the horizon. The frequency of these oscillations is 
more pronounced for the Ghost solution when deeper inside the horizon. In the general 
solution the frequency varies according to the "relative weight" of the Ghost component, 
/So, and the DBI-like one, ao- 



3.2. Power spectrum 

We now turn to the expression of the power spectrum 



2vr2' " 167r(M|,eif2 + 2M|)(ao + v^)='/2|r(| + ^^^^^)|2- 

(23) 

Clearly there is no time dependence in the above result, the modes freeze outside the 
horizon. Note here that, to reproduce standard results, we should re-introduce the speed 
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of sound Cg which, in the language we are using, is related to M2 via: 



1 2M^ 
— 1 - ^ 



cj MlH 



Upon switching to the gauge invariant quantity related to tt by C = —Hn [21] (see also 
the appendix of [36] for the relation at second order) , and reintroducing the proper units 
with Planck mass we get: 

P = (c^o + ^oY'^'H^ .24. 

Again, one could easily check that the above result analytically covers the power 
spectrum of DBI-like and Ghost Inflation. 

As we mentioned in the last comments to equation (14), even when including subleading 
terms in the quadratic action, the functional dependence of our wavefunction does not 
change, only the definition of Q!o,/3o does. Since we are now going to set bounds on 
operators coefficients , we want to be as precise as possible and will therefore extend the 
definition of the two parameters to cover the subleading terms as well. We now have: 

° -Mlfi + 2M| - ?>MlH ' -Mlfi + 2M| - 3Mf//' 

In obtaining Eq. (25), we took into account the fact that the Mi, M2, Ma-driven terms 
multiply operators of the type tt^, (9j7r)^ as well. We also choose to replace M2 and M3 
with a hnear combination of the two masses: we set M| = — 3M| and Mq = M| + M|, 
see also [22]. This procedure allows one to put to zero all the subleading operators 
tuned by M2^M^ and makes the correspondence between inflationary models and the 
switching of the M, M parameters absolutely sharp 

An immediate simplification is that now, upon requiring Mq = <(=^ /3o = 0, one goes into 
DBI inflation, exactly. Similarly, now de-Sitter limit and Mi = give Ghost Inflation 
with ao = 0. The power spectrum looks very similar to the one in Eq. (24) 



4 



p ^ (^0 + ^r'/'H 

^ 167r(M2^eif2 + 2M|-3Mfi7)|r(| + ^^;^^^)r ^ > 



only the definitions of the «0) parameters and the normalization constant have 
slightly changed. Expressing all the parameters in the spectrum in terms of the M, M 
coefficients, we count the degrees of freedom as being five, associated to iJ, e, M2, Mqi ^i- 
The first three are the same that appear also in the standard case as if, e, c^; using Cg 
instead of M2 is just a matter of dictionary. In the Ghost Inflation case, the quantity 
Mo replaces the speed of sound and we are again back to three parameters. In the 
most general case one has to keep both Mq and Mi as well. Bounds can be put on the 
values of these five parameters by employing the expected value for the power spectrum, 
P^. ~ 10"^'^ and its tilt *. Let us also mention here that further mild inequahties must 

+ The reader might worry that one degree of freedom is lost. However, the two coefficients multiply 
basically the same interaction terms in the action up to fourth order. 
* See also [57, 58] 
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be satisfied by the q;o,/3o parameters in order to keep the generahzed speed of sound 
small, Section 2. In Sec. 2.3 below we present a calculation for the tilt and running 
of the power spectrum. These quantities are essentially obtained from considering the 
time dependence of the Mn coefficients in Eq. (26). On the other hand, as has been 
specified above, part of the procedure that led to Eq. (26) has been to disregard the time 
dependence of said coefficients ^. Restoring it at a later step, as we do below, is standard 
accepted procedure because generally only at this stage the effect of time-dependence 
becomes important. A calculation that does without this assumption is presented in 
Sec. 3. 



3.3. Tilt and running 

Below we employ some simplification in order to present our result for the tilt of the 
power spectrum in a way that resembles as closely as possible the typical expression for 
ng — l. Indeed, the spectrum dependence on the Euler F function in Eq. (26) is not to be 
found in e.g. DBI, Ghost Inflation etc. For simplicity, we choose not to write here the 
explicit dependence of the Euler function on the Mn coefficients and leave it implicit; we 
report the full explicit dependence in the Appendix A. We apply the following formula. 



d , „ f dink 
Ug-l^ ——r In Pk = 



dink \ dt 



' 1 dPk 



1 dPk 

t=t* HPk dt 



t=t* J Pk dt 

on the power spectrum, where t* is the time at horizon crossing. The time dependence 
of the M, M coefficients is taken into account and the time dependence of the Euler 
function is dealt with as one would do with a generic function T{t). One obtains: 

f IH^Mle + 8 (2M2^ + 3Mi^) 

~Hf ~ ^ ^ 2 {Hmj,e + 2M24 + SMi^) ^ 
^ 8H^Mj>e'^ - lQH'^Mj>e (2M2^ + 3Mi^) + H'^M],e {-UH^M'^e + 8M2* + 15Mi^) 

2 {H^Mle + 2M24 + 3Mi4) (2H^Mj>e + Mi^ + Mq"^ y^Mm^^TlM^TOM^^ 

{6H^M^e + 24M2^ + 33Mi^) H^Ml>e 



-77 X 



4 {mule + 2M24 + 3Mi4) {2H'^Mle + Mi^ + Mo'^-y/2H^Mj,e + AMi"^ + 6Mi^^ 



H'^Mle Mi 2M2^ 

~^ ^ 4 {H^Mj,€ + 2M24 + 3Mi4) ~ HM2 ^ ]Hm^€ + 2M^^T3M^ ^ 



M2 (6//^M|e + 3Mi^) 2M2^ 



Ml (+3Mi^ - AH^Mle + 4M2^) 3Mi^ 

X 



HMi (H^Mj>e + 2M24 + 3Mi4) (2H^Mj>e + Mi^ + Mq V2^^M|.e + 4M2^ + 6Mi^^ 

ti This is because, at the level of the action, considering the time dependence of these coefficients would 
automatically translate into going at ncxt-to-lcading order in slow-roll; for a leading-order calculation 
it is therefore sufficient to consider their values at the horizon. 
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2 



Mo 6M^ 

X 



HM< 



(2Mo^ + (2// + M,4) ^/=^=£=;) 



Ml SMi^ 



HMi {H'^Mle + 2M24 + SMi^) ' 

where all the quantities are to be intended as calculated at horizon crossing. 
For convenience, we have factored out the usual parameters: e, 77, s (the latter is written 
in M2 language; we give the dictionary in Eq. (30) below) and their generalization: 

Mo Ml f 

''-HMo' ''-HM['' '^-W ^^^^ 

The variable is the one defined in Sec. 2.1 and 77 = e/{He). For simplicity, we have 
defined the variables = M^H^ and = MfH. 

Let us briefly comment on the above results. First note that, as expected, the factor that 
each of these generalized slow-roll parameter multiplies is of order unity or smaller. A 
quick consistency check consists in specializing the formula above to known inflationary 
models, for example requiring Mi = = Mq gives back the usual result [36] for DBI-hke 
models: 

Us — 1 — — 2e — T] — s (31) 

Consider now the more general case with (Mi = , if 7^ , Mo 7^ ), which comprises 
DBI-hke theories and Ghost Inflation models as limiting cases. In such a scenario there 
are some mild bounds to be required on Mq , MpeH"^. First of all, since we are interested 
in the <^ 1 , <^ 1 region of the parameters space this requires: 

M^eH^ < ; < . (32) 

On the other hand, the Mo-driven slow-roll parameters is eo = Mq/[HMq) and therefore 
the above inequalities do not put upper bounds on this slow-roll parameter. Much 
like we will see in the next section for the running of the bispectrum amplitude, one 
must instead be careful to account for the fact that too large a value for Mo could 
give a contribution to the bispectrum amplitude that must be excluded. Indeed, once 
expanded, the Mo-proportional contribution to the quadratic action for the scalar reads: 

J L Q, A J L a d 

(33) 

with the second term on the RHS action clearly contributing to the three-point function 
of the scalar. This Afo-gcncrated contribution must be weighted against the third order 
interactions generated by Mq itself, but also by new Mn coefficients that first appear at 
third order in the action (see [22] and the analysis in Sec. 4)- similar analysis applies 
for ei if we let Mi 7^ 0. 
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We can conclude that the mild bounds on these new, generalized slow-roll parameters 
come from the obvious fact that we are doing a slow-roll expansion, from the value of 
the power spectrum itself and from the requirement that they are not so large as to 

produce too large a value for /atl. 

In order to obtain the running of the power spectrum we proceed as below: 

dus 1 dus , . 

" dh^ ~ H~dt- ^ ^ 

We give below a compact results: 

^ _ ^ I 2 ^0 ^202 £202 vQ-n 



mv^^ H H H H H 

_ _ ^101 _ ^101 _ £000 _ £o0o 
H H H H H ' 

(35) 

and point the reader to the Appendix A for a more explicit expression of the coefficients 
functions ©, ©2, 0??, ©i, ©o- 

3.4- Next-to-leading order 

The discussion presented so far is based on a generalized slow-roll approximation at 
leading order. In particular, the M, M coefficients driving the various operators in the 
Lagrangian are assumed to be time independent. This assumption propagates into the 
equation of motion for the scalar vr, the classical solution itself and the power spectrum. 
To make up for this approximation when calculating the the tilt of the spectrum, one 
restores the time dependence of the coefficients at the level of the power spectrum. 
A more systematic approach consists in accounting for the time dependence of the 
coefficients already at the Lagrangian level, this is done by taking the generalized slow- 
roll approximation to next order. Schematically one has: 

(X J d'^x^f^ - Ml{ZH^{t + tt) + H{t + tt)) + M{t + tt) x (quadratic) . (36) 

The first term does not appear in the Lagrangian in Eq. (14) because, at leading order in 
slow roll, it is not quadratic in fluctuations, it does in fact contribute to the background. 
We see that accounting for the "tt" in the time dependence of M in the last term in 
Eq. (36) would result in a third order operator. This must be considered when studying 
interactions but it is not what we want to analyze here, the wavefunction comes from 
the quadratic Lagrangian. On the other hand, the "tt" in the first term of the action has 
to be accounted for; doing so results in just one additional contribution to the action 
and it turns out to be proportional to e^. This is all consistent with the fact that, at 
leading order, the action is instead proportional to e. We now have: 



M'pH{t){d,7:Y + 2M|(t)7r^ - Mi{t)H{t){2>e - '^^) 
+ ^\{dhr-^MlH{tM. (37) 
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We proceed to write down the equation of motion as in the leading order case, obtaining: 



a, + ao(l + s^)k a, + ^ + s,)^,a, = ^- + Sea H ,^4_^2^_^3ff ) ' 

(38) 



where the following definitions have been employed: 



< 1 ; = < 1 ; f = a\-M'pH + 2M| - SMfH) ; tt 



(39) 

In order to solve the equation of motion one needs to calculate /"// explicitly; the result 
is given in the Appendix A. A compact expression for the equation of motion is given 
by: 

// ~ 2 

+ aok^Uk + Pok^r'^ak = ;^(1 + xo)(7fc , (40) 

where xq is a linear combination of slow roll parameters and q.q,(5q represent a slight 
redefinition of the initial parameters. For explicit expression we refer once again the 
reader to the Appendix A. Equipped with Eq. (40), one uses the Bunch-Davies vacuum 
condition to write down the solution: 

-iaoV^+2^o+^oV9+85^ 1 (2 + ^/9 + 8^^) , -tJJ^k'^x'' 



Q 

a{x,k)^Ci{k) — 



4/3o ' 2 



(41) 

2-i(2+V9+8S^)g-iiV^fe23.2 ^^2-)i(2+y9q^) ^ 

where Q is the usual hypergeometric function, x = — r, and the wavefunction must be 
expanded to first order in Xq = 0. In the first section, we gave the exact expression 
for the solution above at leading order and went on to calculate the resulting power 
spectrum, its tilt and running. The calculation at next-to-leading order has already 
been performed for DBI-like theories of inflation, using the same formalism employed 
here, in [36]. Here instead, we choose to calculate the next-to- leading order Ghost 
Inflation solution obtaining also the tilt of the power spectrum and the running. The 
procedure is a standard one, so we briefly sketch it. The two /c-dependent constant of 
the Ghost equation of motion 

al+pk'r'ak^\{l + x^)ak, (42) 



T 



are reduced to one by imposing the correct leading order limit on the wavefunction. 
The remaining constant is obtained by requiring the proper normalization, that is by 
imposing the following commutation relations to hold: 

[7r(f ), P{y)] = t S^x - y) ; [a,^, a^] = {2nf5\k - p) , (43) 

where P is the momentum conjugate of the scalar tt and the creation and annihilation 
operators a, are the usual operators in the free fleld expansion for the quantized fleld 
TT. Proceeding as prescribed above, one obtains: 

Mk,T) = ^ ^ ^ — . (44) 

3 (1 + f - xo(-| + 7/3 - f + log2)) (r [1 + + r [f] 
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The wavefunction above is the correct Ghost Inflation wavefunction up to next-to- 
leading order in generahzed slow-roll parameters; it gives back the leading order solution 
and respects the proper Bunch-Davies vacuum requirement. From Eq. (44) one can 
readily calculate the tilt of the spectrum: one simply considers the leading behaviour of 
the wavefunction as k goes to zero. In fact, upon expanding for small k one finds that 
our solution goes like: 



4 + y'2^l' + yj^ 



2^2 



k-^-^ , (45) 



from which we obtain that 

- 1 = -4/3x0. (46) 

When specialized to Ghost Inflation, the value of the Xq parameter (which always 
constists of a linear combination of the generalized slow roll parameters) is given by: 

and the running amounts to simply 



.In* =-«°/'^- 

Let us briefly recall what has been done so far at second order in perturbations: we 
obtained a solution to the equation of motion of the tt effective action. This, properly 
normalized, has been used to obtain the power spectrum of curvature perturbations, its 
tilt and running. For the specific case of Ghost infiation we have obtained all of the 
above also at next-to-leading order in generalized slow-roll approximation. 
We now move on to the study of cosmological observables relevant to non-Gaussianity: 
the bispectrum, its running, and the trispectrum. For the bispectrum calculations one 
starts from the cubic effective action as written in Eq. (13). We note here that, for the 
bispectrum calculations, we employed both the exact, general wavefunction as given in 
Eq.(22) and the simphfied solution which is found in models such as DBI infiation or 
K-infiation. By comparing the shape-functions obtained in these two different ways, 
we verified that the simplified solution is indeed a good approximation to the exact 
wavefunction and used only the former in the calculations for the trispectrum. 
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4. Amplitude of the primordial non-Gaussianity: Bispectrum 

In this section we wish to perform a general analysis of the amplitude of the bispectra 
stemming from the general third-order interaction terms. The shape analysis will be 
done in the following section. In the calculations that follow we employ the so-called 
in-in formahsm [26, 28, 27, 29]. To compute the amplitude of the non-Gaussianity, 
indicated by /nl, we proceed as traditionally done in the literature and evaluate the 
three-point correlator in the so-called equilateral configuration where all momenta are 
taken to be equal: ki — k2 = k^. In other words one can write the bispectra of the 
gauge-invariant curvature perturbation ( generated by each interaction term (I) as 

(C(ki)C(k2)C(k3))7 = {27TfS\k^ + k2 + ks)Bi{ki, k2, ks) , (49) 

with an amplitude /j^l defined so that with all three-momenta equal /j^l — 
{6/5)Bi{k,k,k)/P(^{ky, where Pc{k) is the power spectrum of the curvature 
perturbation. For large non-Gaussianitics the linear relation ( = —Htt suffices for 
the bispectrum calculations since quadratic corrections give a negligible contributions. 
Notice also that the values of the integrals appearing in the in-in computations have been 
specified at horizon crossing and the contribution of the integral function at r = — cxd 
has been put to zero mimicking the effect of the slight rotation of the r-axis into the 
imaginary plane. 

Given the broad number of possibilities, we choose to compute numerically the 
amplitude of the various bispectra identifying six benchmark points in the {ao, Pq) plane 
and numerically integrating the exact wavef unctions. Since one is interested in probing 
models with large non-Gaussianity, the values of (q;o,/3o) are taken much smaller than 
unity as described in the following Table: 



Benchmarks 


1 


2 


3 


4 


5 


6 




10-2 





0.5 • 10-2 


2 • 10-^ 


10-^ 


10-6 







0.5 • 10-^ 


0.25 • 10-^ 


5 • 10-5 









Upon using Eq. (18) one can check that the first four benchmark points correspond 
to the same choice of the effective horizon. This choice has been made to suitably 
perform a comparison between the various cases and in particular against the values of 
/nl of purely P{X, 0) and Ghost models which correspond to benchmarks 1 and 2. The 
last two configurations (which are P{X, 0)-typc) probe the space of extremely small ao 
and, in interactions with at least two space derivatives, are expected to give a larger 
amplitude than the first four points, at least if the interaction terms are regulated by 
unconstrained masses. To get the feehng of the figures involved, if we restrict ourselves 
to the case of theories for which ao reduces to the usual sound speed c^, typical values 
of Cg are between 10"^ and lO"^. Note also that the definition of varies according to 
which operators are switched on in the action (one can well be in the de Sitter limit 
where the only spatial quadratic term has four derivatives; this leads to a different c^. 
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see also [31]). Therefore we choose to specify all the amphtudes as a function of the 
various M and M masses. 

It is convenient at this stage to clearly point out which M coefficients in Eq. (13) 
are free and underline the relations among the constrained ones. From Eq. (25) one can 
see that, despite fixing ao, as we did in the Table, there is still the freedom to pick any 
reasonable value for either M| or Mf. Similarly, fixing Pq does not completely specify 
Mo- In other words, both Mf and Mq are constrained by our choice of the (etc/So) 
parameters; all the other coefficients are unconstrained. As the M's are expected to 
set the energy scale of the various underlying theories, they should be larger than the 
Hubble rate H, and can go up to Mp\. As elucidated in Ref. [49] though , for the action 
(13) to be as general as possible, one might want to require the M's to be smaller than 
the Planck mass. That said, some useful inequalities that the mass coefficients must 
respect can now be reminded. Due to the fact we are working in the decoupling regime, 
we must require M| < Mp^H^. Also, the fact that we are probing the {ao, (3q) <^ 1 space, 
imposes bounds on some masses. Consider the parameter ckq in Eq. (25). There are two 
ways this coefficient can be much smaller than unity. The first and perhaps most natural 
way, is to ask M| » Max (— Mpjij, —MfH) which, due to decouphng inequalities on M2 
puts a bound on Mi, Mf <^ Mp^H. The other possibihty requires a partial cancellation 
in the numerator of ao, —Mp^H ~ MfH/2 which is certainly possible but it implies we 
are neither in the DBl, nor in the ghost regime, both of which have Mi = 0. Looking 
at (3q we see it is enough to require MqH'^ <^ or MqH"^ ^ —MfH and again, the 
first condition seems more natural. Let us stress here that, upon requiring the masses 
to be all of the same order, M, and using that H/M <^ 1 in the effective theory, one can 
easily obtain small ao,(3o coefficient. However, when employing a single mass scale M 
in the whole Lagrangian, working with tiny values for cto and ^0 would put a bound on 
M and necessarily influence the magnitude of all the interaction terms. In our analysis 
we let the M, M's coefficients be not all of the same order (with some important caveats 
upon which we expand at the end of this section). 

Below we present the results for each interaction term. All the amplitudes can 
be written as a dimensionless coefficient, 7„, times an {ao, /3o)-dependent numerical 
coefficient. The terms described in the first subsection are interactions that have already 
been discussed in the literature. The novelty here is represented by the fact we are able to 
study also interpolating configurations through the third and fourth benchmark points. 
In the second subsection we report on the amplitudes of the contributions from some 
curvature-generated terms that have never been discussed in the literature. 

4.I. Amplitudes from P{X,(f))-type interactions and first two curvature- generated terms 

The amphtudes from "DBI-hke" (an expression which we use as a synonym of P{X, 0) 
models) interactions and first two curvature-generated terms are the following: 

• Oi = -2Mln{dinf/a^ 
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benchmarks 


1 


2 


3 


4 


5 


6 




10^71 


8 • 10^71 


6 • 10^1 


4 • 10^71 


10^71 


10^71 



where 



7i = 



Ml 



2M| + M^^eH^ - 3MfH' 



(50) 



We see the parameter 71 can in principle be of order unity. Indeed, if one assumes M| 
is the largest term in the denominator (in DBI this would correspond to a very small 
speed of sound), 71 is roughly 1/2. 



O2 = -4/3 M^ir^ 



benchmarks 


1 


2 


3 


4 


5 


6 


fMa 

Jnl 


1/272 


10-272 


5 ■ 10-272 


10-272 


1/272 


1/272 



where 



+ M^^eH^ -3MfH' 

The parameter 72 can be even larger than unity if, for instance, M3 is larger than M2. 
We can see, already at this stage, the effect of small values of and /3o at work: the 
numerical factor of a spatial derivative-free interaction is much smaller than that of a 
third order term like the M2 one calculated above which has two spatial derivatives. 



O3 = -1/2 MUdiTrfdH/a' 



beiichiiiarks 


1 


2 


3 


1 


5 


6 


fMi 

Jnl 


10^73 


10^73 


4 • 10^73 


1.5 • 10^73 


10^73 


lO^S 



Due to the fact that Mf is a constrained parameter we find that the amplitude equals 
an (cco, /5o)-dependent number times the paramter 

M^.eH^ + 2M| - SMfH ' ^ ^ 

One should substitute the various values of cto in the Table in the 73 expression given 
above. Barring cancellation between different mass terms, 73 is generally smaller than 
one and can be as small as ckq itself. 



04 = -1/2 MoV4 i5H{df7:){dj7:y + ijrdfdjTrdjT:) /a' 



benchmarks 


1 


2 


3 


4 


5 


6 


Jnl 


3 • 10^ 74 


10^74 


1.3-10^74 


10^74 


3 • 10^ 74 


3 • 10^3 
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with 

74 = {MlH'^)/{Ml^eH'^ + 2M^ - SM^H). (53) 
We see that large numerical coefficients appear in such a case, nevertheless 74 = /3o <^ 1- 

4.2. Amplitudes from curvature- generated novel interaction terms 

We come to the curvature-generated interaction terms that generate novel bispectra. 
Their amplitudes are given by 

• C»5 = -2/3 Miir'^dfTr /a^ 



benchmarks 


1 


2 


3 


4 


5 


6 


/•M4 
/nl 


7-10275 


10^75 


2-10275 


10275 


7-10^75 


7-10^75 



where 

75 = (M^H) /(M^iei/2 + 2M| - SMfH) . (54) 

The coefficient 75 can be larger than unity. Barring cancellations in the denominator, a 
75 3> 1 larger than unity imposes M4 ^ M2, Mi. 

• e»6 = 1/3 MMdfnf /a^ 



benchmarks 


1 


2 


3 


4 


5 


6 


JNL 


5-10^76 


1.6-10^76 


2-10^76 


1.6-10^76 


4-10^76 


4 - 10^2^^ 



where 

76 = {M^H^)/{M^^eH'^ + 2M^ - SMfH). (55) 
To get 76 S> 1, one needs to impose a less natural condition M^H^ 3> M|, MfH. 

• O7 = 1/3 Min{dy7if /a^ 



benchmarks 


1 


2 


3 


4 


5 


6 


/•Me 


10^77 


4-10=^77 


5 - 10^ 77 


4 - 10=^ 77 


10^77 


10^2^, 



where 

77 = (M|i/2)/(M2jei/2 + 2M^ - 3M^H). (56) 
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The same consideration as for the case of M5 apply here. Note again that the numerical 
values, especially in the fifth and sixth benchmark points tend to be much larger for the 
M coefficients with the most spatial derivatives, thus conferming our expectations. 



• C>8 = -1/6 Mrid^irf /a^ 



benchmarks 


1 


2 


3 


4 


5 


6 


JNL 


8-10678 


2.6 • 106 78 


3.5 • 106 78 


2.6-10678 


8 • 10^278 


8 • 10^^78 



with 

78 = (M7H^) /{M^ieH^ + 2M^ - 3M^H) . (57) 

Notice that in this case , and in the following ones, the 7 coefficients are naturally ex- 
pected to be smaller than unity. For this interaction, as well as the two following ones, 
the numerical factor coming from the integration can be quite large, especially in the 
fifth and sixth benchmark points. This is clearly due to the six space derivatives that 
characterize these interactions. 



benchmarks 


1 


2 


3 


4 


5 


6 


JNL 


2 • 10679 


6 • 10^79 


8 • 10^79 


6 • 10^79 


2 • 10^279 


2 • 10i%9 



with 

79 = (MsH^) /{M^i^H^ + 2M2^ - 3M^H) . (58) 
Let us stress that one could guess the amplitude of M^^g terms by simply looking at the 
results obtained for the interaction term tuned by M7 because, althought these terms 
might produce a different shape for non-Gaussianities, they have essentially the same 
structure as far as the integration is concerned. 



• Oio = -1/6 Mg dijirdjkTrdkin jd 



benchmarks 


1 


2 


3 


4 


5 


6 




106710 


3 • 10^710 


4 • 10^710 


3 • 10^710 


10^^10 


10^^10 



where 

710 = {MgH^)/{M^^eH^ + 2M| - 3Mi^i^). (59) 

One could again just read off the maximum value from /j^^, they differ by just a factor 
1/8. 
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4-3. Some general considerations on the amplitudes of non-Gaussianities 

It might be worth now to pause and comment on our findings in relation to earher results 
in the literature. As we have stressed, the non-Gaussianity generated from operators 
proportional to the M's masses arises from curvature terms that are often neglected. 
From general considerations, one expects that the curvature terms, especially the ones 
coming from SK^-type of contributions, can often be neglected. This is not generically 
true though. The general structure of the amplitudes we have found for the bispectra 
can be schematically written as 



where Mn stands for a generic mass coefficient of the type M, M. 

If Mn ^ H/(ao + (c^o + S/^o)"*"^^)"*"^^, and the masses are all of the same order, then 

it is natural to expect the Mn with the largest exponent to dominate. In this case 

the dominant terms are those associated with M2 and M3, corresponding to DBI-like 

infiation. Furthemore, we may also recall that large masses appear also in the definition 

of ao, Po and, at least in P{X, 0) models, a large M2 is needed to have a small speed of 

sound. 

On the other hand, one can employ the freedom for all the M's not to be of the 
same order in magnitude. The Mn's have a natural upper bound that must be smaller 
than the Planck mass and in the general theory with both DBI and ghost (and more in 
general curvature terms) operators switched on, one can allow for a very small speed of 
sound without assuming much on the unconstrained masses. In such cases, the general- 
ized sound speed can be so small that H/{ao + (ckq + 8^0)^''^)^''^ is actually larger than 
Mn. Consequently, the amplitudes of the masses with a smaller exponent, typically the 
barred M's, are not negligible any longer. Wc conclude that large non-Gaussianities 
may be induced by theories parametrized by suitable values of the M masses. 

It is important at this stage to offer some additional comments on the general 
structure of the Lagrangian and of the resulting amplitudes. As clear from our discussion 
in Section 2, each space derivative acting on the scalar tt can be schematically written 
in Fourier space as a term Hir divided by a number which is much smaller than 
unity in the cases of interest. Now, it is clear from dimensional analysis that other 
terms in the action, specifically the ones with less derivatives, will have their own 
if/(ao + (cKo+S/^o)^^^)^''^ factors coming from derivatives but, most importantly, will also 
need some generic M factor to make the action properly dimcnsionless. This is exactly 
what happens in our theory (sec Eq. (13)) as well: the exponent of Mi is bigger than the 
one of M2 and so on. Wc have employed in the analysis of the amplitudes some freedom 
on these M„ parameters to show that even terms with higher spatial derivatives can give 
non negligible contributions to the amplitude. Of course, if one wants to have a reliable 
effective theory the Mn coefficients must eventually prevail over H/{aQ + {al+S/3oy/'^y/'^ 




(60) 
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(in a Lorentz invariant theory it would suffice to ask for M„ ^ H, here we need more) 
so that it makes sense to consider higher derivatives up to some given finite order, but 
not further. 
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5. The shapes of non-Gaussianities: Bispectrum 



In this section we wish to analyze the shapes of the bispectra generated by the various 
operators analyzed above. In calculating the amplitudes of the bispectra in the previous 
section we have chosen to perform, albeit numerically, calculations with the exact 
wavef unctions. We decided to counterbalance the loss of information in not having 
ao, explicit in the result by running the same procedure for six different benchmark 
values of the parameters. In analaysing the shapes we do not enjoy the possibility 
to use the exact wavefunction any longer as we must approximate the wavefunctions 
inside the conformal time r integral(s) calculated as prescribed by the in-in formalism. 
Of course, exact results arc always available when the classical solution reduces to the 
usual, Hankel function H(2,/2){t) ■ In computing terms like 



we have choosen to expand in series each /c-mode inside the integrand within a region 
that starts from slightly inside its effective horizon (setting to zero the function in the 
rest of the interval). This enables us to keep both, parameters and external momenta, 
arbitrary. This approximation is justified by the fact that, due to the oscillatory 
behaviour of the wavefunctionsff, the main contribution to the integral comes only 
from the region where all the wavefunctions are not oscillating anymore (to be safe, we 
actually choose to include the region where the /c-mode with the latest effective horizon 
is still within its horizon). Furthermore, it is reasonable to assume that, even if some 
non-negligible contribution is being left over because of this approximation procedure, 
it might have a systematic effect on all A;-modes and would therefore not change the 
shapes of the bispectrum of perturbations. This procedure has been calibrated with the 
calculations of Ref. [51] which were performed exactly: we are able to reproduce the 
very same shapes for the bispectrum. 

Following the scheme used for the amplitudes, we now report below the shape of 
non-Gaussianity for each interaction term. Within a single interaction term, we consider 
four configurations in the (ccq, /5o) -plane. We preliminary found that the shape are not 
particularly sensible to their absolute values. This is expected as these coefficients, 
being the same for each fc-modc, enter mainly in the amplitudes, not in the shapes. 
Therefore, we limited our attention to the Q;o//3o ratios. We probed the following 
cases: cto 7^ 0,/3o — (DBI-like); — 0,/3o 7^ (ghost); ^ Q ^ (5q (with the 
solution interpolating between DBI-like models and ghost inflation and «o playing the 
dominating role in determining where is the effective horizon, we call this configuration 

ft Of course, an oscillatory behaviour is not, by itself, enough to provide a cancellation all over the 
-inside the horizon- region. Indeed in some cases for r — >■ —00 the amplitude of the wavefunctions 
increases and one certainly docs not expect this to give zero contribution. On the other hand, much 
like in the simplest single field slow-roll calculations, one is expected to slightly rotate r to the imaginary 
plane to match the vacuum thus basically putting to zero the contribution at —00. 




(61) 
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^4) and 7^ 7^ /3o (with the effect of f3o leading in the expression for the horizon, we 
call this configuration B). 

Let us briefly add some general preliminary comments. We need to point out 
that when we use the name DBl might generate some confusion. Indeed, we use it 
also to describe shapes due to interaction terms coming from curvature perturbations, 
such as Ml 2,3..... What we mean here is that we can employ in the integrations the 
usual wavefunction Tr^ oc e~*v^^'^(i + iy^A;r) without resorting to any approximation. 
This wavefunction can be used as far as it is the solution to the equation of motion. 
This is a condition concerning only second order perturbations: we need only require 
/3o = 4^ = Mq to employ it. Wc can, at the same time, have extrinsic curvature- 
driven interaction terms, the M's, and yet find exact results. The shapes corresponding 
to what we denominate DBl configuration have been thoroughly investigated in a 
number of papers [18], but what has been generally left out is the contribution coming 
from the extrinsic curvature terms: it is indeed possible to have the usual Hankel, 
-f^3/2(T), wavefunction as a solution to the classical equations and, at the same time, 
switch on curvature operators like Misgjgg. All the shapes obtained in the DEI 
configuration are generated through exact analytical methods. We will call them 
as "exact-DBI configurations". Should we find, as we will, a shape which is not 
equilateral in the first configuration, that shape suffers none of any possible limitations 
the approximated method might introduce. The ghost configuration has been analyzed 
in depth in many articles, among which [35, 31], and again, not all the terms coming from 
extrinsic curvature have been taken into account. It is important to note though that in 
[21, 35, 31] the M coefficients multiplying curvature terms have all been chosen of the 
same order thus resulting in only a couple of leading curvature terms (they correspond 
to the Ml, Mo contributions). Finally, the shapes in configuration A, B have never been 
analyzed before. 

The shapes reported below have been obtained by employing the shape function 
B(1,X2,X2), where Xi = ki/ki, which posseses the same fci 2,3-dependence as the three 
point function. What is plotted exactly is X2xfB{l,X2,X3)/B{l,l,l) in the region 
satisfying X2 > x^ > 1 — X2. Let us remind the reader that a shape is called local 
when it peaks for a small value of, say, X2, with Xi = X3 ~ 1; equilateral when it peaks 
in the equilateral configuration xi = l,a;2 = 1 = 2:3 and is called flat for squashed 
triangles with Xi — 1 — X2 + x^. In particular, we find, as detailed below, that some 
novel curvature-generated terms produce a flat bispectrum which specifically peaks for 
3^2 = 3^3 = 1/2 = l/2xi. In presenting the shapes we follow the same order and 
organization we employed for the amplitudes: 

5.1. Shapes from P{X,(f))-type interactions and first two curvature- generated terms. 

• Oi = -2M|7r(9,7r)Va2 

Note that, although more or less sharply, all the four plots are peaked in the equi- 
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0.0 



Figure 3. DBI configuration on the left, obtained using exact methods; the 
approximated ghost shape on the right. 




0.0 0.0 



Figure 4. A configuration on the left, B configuration on the right. 



lateral configuration. 
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C»2 =4/3 M|7r3 





We obtain equilateral shapes in all four cases. It is somewhat expected that the 
general, interpolating solution employed in configurations A and B, will give qualita- 
tively the same plot, we have verified it in these first two rounds of shapes. 
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.03 = -1/2 Ml{diTifd]n/a^ 




Figure 7. exact DBI configuration on tlie left; approximated ghost shape on the right. 




0.0 0.0 



Figure 8. A on the left, B configuration on the right for the Mi-driven interaction 
term. 

Also the bispectrum generated by this interaction term has an equilateral shape; 
the last two plots show, employing the general wavefunction, that also the interpolating 
models produce an equilateral shape. 

. C4 = -1/2 MoV4 (5iJ(927r)(9,7r)2 + AiidjdjT^djTi) ja^ 

The various bispectra peak in the equilateral configuration. 
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1.0 




Figure 10. A on the left, B configuration on the right for the AIq interaction term. 



5.2. Shapes from curvature-generated novel interaction terms. 



• Cs = -2/3 Mi-k^dfTc lo? 

With the term tuned by M4 we start including in our description the contributions 
that have so far been neglected in the literature. 





Figure 11. exact DBI configuration on the left; approximated ghost shape on the 
right. 
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Figure 12. A on the left, B configuration on the right for the M4-driven interaction 
term. 

Equilateral shapes for all configurations are obtained. 



Oe = 1/3 Ml7r{df7cy /a^ 





Figure 13. exact DBI configuration on the left; approximated ghost shape on the 
right. 





Figure 14. A on the left, B configuration on the right for the Ms-driven interaction 
term. 

Again, all the four bispectra peak in the equilateral configuration. 
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• 07 = 1/3 Mlnidi^nf /a^ 




Figure 15. exact DBI configuration on the left; approximated gfiost sliape on tlie 
right. 




Figure 16. A on the left, B configuration on the right for the Mg-driven interaction 
term. 

This is one of the interesting novel curvature-generated terms that give rise to a 
flat shape (more precisely, the plot peaks at fci = 1, ^2 ~ 1/2 ~ k^)- Note that for a 
very similar interaction, namely the one generated by Oq, we saw an equilateral plot. 
Here, derivatives combine to provide a different /c-dependent factor outside the integral. 
Writing in Fourier space the interaction term Oq we obtain something proportional to 
fc2^^3^, while here we obtain a contribution proportional to {k2 ■ ^s)^. 
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= -1/6 Mridfnf ja^ 





Figure 17. exact DBI configuration on the left; approximated ghost shape on the 
right. 





Figure 18. A on the left, B configuration on the right for the My-driven interaction 
term. 



All shapes peak in the equilateral configuration. 
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Figure 19. exact DBI configuration on the left; approximated ghost shape on the 
right. 




1.0 x3 1.0 



Figure 20. A on the left, B configuration on the right for the Mg-driven interaction 
term 

This is a second interaction term that produces, just as for O-j, a flat shape for the 
bispectra. Comparing it in Fourier space with our findings for Og? one can see that it 
is due to the way the spatial derivatives are combined. As shown in Figs. 21 and 22, 
also Cio gives rise to flat-shape bispectrum. We can see that the interactions Cs, 
and Oio have the same structure as far as the integral is concerned; on the other hand 
their /c-dependence goes like kl /c| /cf (/c2 ■ k^Y + perm and {ki ■ k2){k2 ■ k^){k^ ■ ki), 
respectively. The last two produce a flat shape. 
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• Olo = -1/6 MgdijTTdjkndkiTT /a! 




Figure 21. exact DBI configuration on the left; approximated ghost shape on the 
right. 




Figure 22. A on the left, B configuration on the right for the Mg-driven interaction 
term. 

The shapes are peaked in the flat configuration. 

Obviously, having presented all the shapes due to each individual interaction, one 
might proceed with the study of the shape of linear combinations of them, much in 
the spirit of the orhtogonal shape recently introduced [31]. On the other hand, we 
are using approximated methods for three of the four configurations under scrutiny 
and it is therefore not a cautious step to infer new shapes from linear combinations 
of approximated ones, especially when delicate substractions are involved. One could 
proceed to study linear combinations in the case of the first configuration. We leave 
this to future work. 

5.3. Some general considerations on the shapes of non-Gaussianities 

From the above discussion we can read off some general qualitative features 
characterizing the shapes of the bispectrum of a generic single-field model of inflation: 

• Two qualitatively very different kinds of shapes appear: equilateral and flat. As 
pointed out before, the next natural step would be to consider the shapes resulting 
from linear combination of the various interaction terms in the spirit of Ref. [31]. 

• In single field-models of inflation usually large non-Gaussianities are associated 
to equilateral shapes. In order to obtain a flat shape one needs to consider linear 
combinations of interaction operators such as what is done in [31] or models with an 
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initial vacuum different from Bunch-Davies. Interestingly, in our case a flat shape 
emerges simply from individual operators generated by curvature-related terms. 
When the fiat shape appears, it does so in all four configurations considered. We 
stress this point because it implies that this result does not depend on the type of 
wavefunction one employs in the calculation, be it the DBI inflation solution, the 
ghost inflationary one and the exact wavefunction that covers also more general 
models. The results of the DBI-like configuration are exact and easily reproducible 
with analytical methods. In fact, the DBI wavefunction is the usual solution of 
the standard single-field slow-roll inflation. Adopting such a wavefunction we may 
provide analytic results for the flat bispectra 



_ {60a%^ + lla%^ + b^ + (48a^ - Aa^b^ - 36^) - Ab^c'^) H^M^ 



_ (15a^ + b^ + 3bc^) {12a^ + b^ - 6b^c^ + 86^ c^ + 8a^ {b^ - 2bc^)) H^Ms 



_ 3 (15a^ + b^ + 3bc^) (8a^ + b^ - Gb'^c^ + 8b^c^ + 8a^ (b^ - 26c^)) H^Mg 

{CklCk2Ck3}M, - 2^X^WW| ' 

(62) 

where the overall momentum conservation delta has been omitted, 

a = (/ci/ca/cs)^/^ b^ki + k2 + k3; c ^ {kik2 + kiks + k2k3y^^ , 

and CKo is now the usual speed of sound: 

ao = -Ml^H/{-Ml^H + 2M|) = . 

Note that all the three expression given above have a maximum precisely in the flat 
configuration [ki — l,k2 — 1/2 — /ca). What immediately stands out in Eq. (62) is 
the presence in the numerator of factors consisting of sTibtractions between generally 
positive /c-symmetrized terms: it is this characteristic that selects a flat, rather 
than a equilteral shape, as one can readily verify by checking the bispectrum of the 
"equilateral" interaction terms. 

The expression in Eq. (62) is exact for all P{X, (p) models and can be employed 
to get a shape qualitatively similar in the other three configurations. For practical 
purposes, we give below a very simple expression that very closely mimics the 
behaviour of the typical bispectrum contribution that generates a fiat shape: 
, , , i-kf + k2^ - k^) {kl + k2^ - k^) {-kl + k2^ + k^) 

^^''^'^^^'^^ klk2^k,^h + k2 + k,y ■ ^^^^ 

As a general rule, the terms which are going to generate a flat shape can be read 
off already at the Lagrangian level: indeed the flatness originates from the way 
the external momenta combine with each other and are summed over. Whenever 
mixed space derivatives act on a single tt term and the mixing is repeated on at 
least another tt field, the shape turns out to be fiat (note that this criterium puts 
M6,8,9 contributions in the same, "fiat" class, but correctly excludes apparently very 
similar ones such as M2,3j). 
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6. Running of f^L 



In the quest for properties that help in removing degeneracies among the many 
inflationary models one generally considers another observable beyond the power 
spectrum and its running, i.e. the analysis of non-Gaussianities. Starting with the 
bispectrum, one can study its amplitude, shape and running. In the same spirit of the 
analysis we performed for the power spectrum, we now want to estimate the value for 
the running of the bispectrum amplitude, Jnl- In performing the calculation for P^, 
Ug — I and then a^, we used the fact that the coefficients driving quadratic operators in 
the Lagrangian are nearly constant, up to slow-roll corrections. At first approximation, 
one writes down the power spectrum as a function of these parameters calculated at 
the horizon. Only when calculating the tilt of the spectrum and its running one does 
consider the time dependence on the Mn's, thus obtaining Eq. (28), (35). Similarly 
here, we will employ the results on the bispectrum amplitude contributions generated 
by independent interaction terms as given in [22]. More precisely, we are going to focus 
on the running of f^L generated by a third-order interaction term whose bispectrum 
shape-function peaks is an uncommon fiat configuration (see also [31]). In fact, as wc 
have seen in the previous section and as detailed in [22], there are several independent 
terms that generate at least two qualitatively different ffat shape-functions; we choose 
here to concentrate on the analysis of just one term as the same considerations and 
conclusions can be straighforwardly adapted to all of them. 

The analysis for the running of Jnl has been done for several inflationary mechanisms 
such as DBI inflation, and others. We are going to extend this type of study to 
a third order interaction term driven by a nearly constant coefficient, Meit), and 
which is generated by an extrinsic curvature contribution which we reproduce here for 
convenience: 

M6(t + 7r)2 



/ 



d^xdt^J—g 



... + 



TT 



+ 



(64) 



We have already calculated the corresponding contribution to f^L in six different points 
of the (cco, /5o) plane, we report them in the Table below. 



Benchmarks 


1 


2 


3 


4 


5 


6 




10-2 





0.5 • 10-2 


2 ■ 10-^ 


10-4 


10-6 







0.5 • 10-4 


0.25 • 10-4 


5 • 10-5 









Configuration 1 describes pure DBI-like theories, pure ghost corresponds to configuration 2. In 3,4 a more general model 
is considered while in the last two configurations one aims at considering the cases characterized by very small generalized 
speed of sound, -^/oo- 



Corresponding to each one of the configurations of the table above, we report below the 
values for /jvl- 
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benchmarks 


1 


2 


3 


4 


5 


6 


Jnl 


10^76 


4 • 10^ 76 


5-103 76 


4 • 103 76 


10^76 


10^2 



The value of the bispectrum amplitude fj^^ in the six different configurations. 



The dimensionless quantity 75, which was introduced in Section 4, is given by 

76 = {M^H^)/{Mj>ieH'^ + 2M^ - 3MfH). (65) 

Computing the running of Jnl amounts then to calculating 

d\n\fNLik)\ 1 dfNL 

dlnfc "i//., dt ■ ^''^ 
For the specific case at hand, we have: 

1 pMg I d [ 



HfSl Hf§ldtl ^ ' ^1 t=t* {H'^Mle + 2M2^ + SMi^) 



4M2^ + 2Mi4 



+ 2e6 + e;v + 



H'^Mle 8M2^ 4Mi^ 



^ii"2^|,e + 2M24 + 3Mi4 (ii'2M|,e + 2M24 + SMi^) (ii'27v^2e + 27^2^ + SMi^) ' 

(67) 

In the above expression N stands for a numerical factor which is dependent on and 
/3o, specifically it goes like (ag + /^o)"" with n small positive integer. The generalized 
slow-roll parameters are defined as usual, with the only new ones being cq^cn'- 

66 = ^^; ejv = 0(1) X (6,77,62, ei,eo). (68) 

The last expression in Eq. (68) means that can be expressed as a linear combination 
of the generalized slow roll parameters introduced before multiplied at most by an order 
unity constant. Overall, we can then conclude that the running of f^^l is of the order of 
the generalized slow roll parameters (or a linear combination thereof). Now, all these 
parameters, except for 66, can be expressed as functions of M, M coefficients that first 
appear in the quadratic Lagrangian of the theory. As such, four of these parameters 
could in principle be expressed as a function of the observables one usually uses, namely 
the scalar and tensor spectral indices and their running. Not so for 65, on which, in 
principle, we enjoy more freedom as it drives terms in the action that are at least 
cubic. It could indeed be that the running of /^^, and of Jml itself, is dominated by 
this contribution and could therefore be larger than one finds in some general single- 
field slow-roll models where the extrinsic curvature-generated interaction terms are not 
accounted for. The same considerations apply to other interaction terms such as the 
Mg-driven one. On the other hand, special care must be exerted so as to make sure that 
requiring 66 to be the leading generalized slow-roll parameter in Eq. (67) does not spoil 
the possibility to have the contribution in Eq. (64) dominate the overall bispectrum 
amplitude, which is what made this type of contribution interesting in the first place. 
Indeed, there are two ways of making the leading parameter, a large Mq and a small 
Mq. Pushing the latter option too far the interesting and possibly leading bispectrum 
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and trispectrum amplitudes [22] would become subdominant and this would make the 
corresponding flat shape-functions a mere curiosity. But also the former option has to 
be discussed and this is clear from the expansion of Eq. (64) resulting from considering 
the Mq time dependence: 



d'xdt,/^[... + J2(d,jnf + J2i^^J7^f x 7r...J . (69) 



Indeed the quartic interaction term that appears is proportional to Mq and if the latter 
is too big it could give rise to too large a contribution to the power spectrum at one 
loop and would have to be ruled out. Below we give a number of inequalities that 
the quantities Mq, Mq need to satisfy in order not to spoil the appealing bispectrum 

and trispectrum features outlined above. It turns out they arc not too restrictive and 
that the running of f^l can be safely ruled by eg. We first write down the inequalities 
that stem from requiring that the Mg-proportional quartic interaction term is not the 
leading interaction in the fourth order Lagrangian (this would have consequences on 
the loop corrections to the power spectrum as well) as compared to the usual single- 
field interactions (for the complete action at fourth-order see Ref. [23] or the detailed 
derivation in the next Section): 

M| (ao + x/^)M| {a, + V^)^Mi 

>eQ:^e, >eQ:^e, >eQ:^e. [{0) 



Notice that only one of these inequalities need be satisfied. The eg » e part of 
the inequalities above ensures that indeed eq is the dominating generalized slow-roll 
parameter. One has to keep in mind here that large non-Gaussinities are generated by 
requiring the generalized speed of sound, ~ ao + \/A)) to be much smaller than unity and 
so the first of these inequalities seems somewhat less stringent than the others though 
there is no requirement on the Mn's to be all of the same order. 

As anticipated a small Mq can in principle lead to a subleading contribution to the 
bispectrum signal thus rendering the corresponding fiat shape- function less interesting. 
Borrowing the third order action of Eq. (13) and employing the estimates on the 
wavef unction we showed to hold at horizon crossing, one is able to derive the inequalities 
below: 

>i; ..A. . /7^x9 >i- (71) 



M|(ao + v^) ' M|(ao + v^)2 
This time it is necessary that both the inequahties in Eq. (71) are satisfied. The 
coefficient M4 is not found in the equation above as it first appears in theory Lagrangian 
at fourth order, so it is not involved in the tree level bispectrum calculations. The further 
freedom on M4 (and other cocfTcicnts) that results from this simple fact can be used to 
study models of inflation which present a relatively small bispectrum together with a 
larger trispectrum signal [23, 59]. The inequahties given in Eq. (70) and (71) above are 
indeed compatible for values of (ao + ^/Po) smaller than unity, an assumption which is 
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generally made when looking for models that can produce large NG, as shown in [22, 23] 
and detailed in Section 2 and 3. We see then that there is a whole, large window of 
values for the coefficient Me that would allow for a running of Jnl dominated by ee- This 
effect is generated purely by third order terms and specifically by extrinsic curvature- 
generated interaction terms driven by coefficients which first appear in the Lagrangian 
at third order. 

Let us consider a specific example with realistic values for the generalized speed of 
sound. For simplicity, we take ccq = 10~^ ; = 0: this corresponds to a generalized 
speed of sound ~ y/a^ of about 1/10. If we want the corresponding value of the Me 
contribution to /tvl to be well within the WMAP7 [15] limits for say, f^^"^"^^^ * , 
then one automatically obtains from the first column of the Table above Eq. (65) the 
inequality 7e < 2 x 10"^. Our first inequality in Eq.(70) is easily satisfied in the case at 
hand and for the first relation in Eq.(71) one obtains 



A/J(a„ + v^) 10-2 

So that, indeed, there is room for a leading ee in the running of Jnl also after the bounds 
on Jnl have been duly set into place. 

Adding the results of this section to the analysis of [22, 23] one can safely say that, with 
respect to all observables one is ultimately interested in, the analysis of these curvature 
terms has shown they can have leading effects on all quantities and must therefore 
always be included in a thorough analysis of non-Gaussianities. 



*Implementing the bounds such as —410 < f^^°^°'^^^ < 6 is indeed the best we can require 
on the Me-driven contribution to fj^L as the shape-function generated by Mq is certainly closer 
to the orthogonal shape than to the equilateral or local one. 



50 



7. The Hamiltonian up to fourth order 



Before we proceed to write down the quartic Lagrangian, let us tress that the cubic 
Lagrangian is going to be just as relevant for the four-point function calculation as it 
contributes to the so called scalar exchange diagram we will shortly calculate below. 
Adding the trispectrum analysis to the bispectrum data we gained in the previous 
section, here the hope is to be able to identify distinctive features for as many as possible 
different combinations of the M^'s in the form of specific patterns they produce in the 
shapes of the various correlators of curvature perturbations. The degeneracies among the 
results for different inflationary mechanisms that will inevitably arise might be removed 
by a joint analysis of the different n-point functions, starting with the bispectrum, the 
trispectrum, loop corrections to the power spectrum and so on. Let us briefly go through 
some of the main features of the third order effective action above. All the comments 
can be straightforwardly extended to the fourth-order expression as well. In deriving 
the fourth-order Lagrangian we use the same algorithm that was used at third order so 
there is no need to reproduce it here. We do a small exeption for the following estimates 
we showed to hold at horizon crossing: 

TT ~ i/TT, Vtt ~ — -j^^^^=^^^^= n = H/cs n, (73) 

because they will turn out to be, once again, very useful. The most general fourth-order 
action in the usual set up is: 



/ 
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+ 



4! a6 



^ Mu{t)27rdl7r{dij7rf ^ M^^{t) 2 iidijTidjkT^dkiT: Ni{t) {dfTr)^ N2{t) {dlnf {di^Tif 



4! 



4! 



4! 



4! 



4! 



4! 



4! 



(74) 



Note that, as pointed out in [32], starting at fourth order in perturbations, one cannot 
immediately read off the Hamiltonian from the expression of the Lagrangian, in other 
words H — —L does not hold here. We use the results one obtains by adopting the 
correct procedure which was outlined in detail in [32]. 

Let us split the interaction Hamiltonian we will be concerned with as — + H^, 
one can prove that the overall interaction Hamiltonian is then: 

1 



int 



d a:v— (? 



2M| + H^eMl - ?,HMl 



2a8 



3a4 



3a^ 



3(f 



3a/ 



47r2(927r)2M| {dinfidlirfMlMi 27r\dfTTfM^Mi {dlnfdfdjTrdjirM^M, 



9a^ 



3a6 



3a^ 



6a^ 



{dlTrfdjTTdjTrMiMi 2i:{dl'KfMlMl {OIt^YMI {dk7rf{dij7rfMlM, 



6a^ 

27r^{dij7rfMiM, 



9a6 



+ 



36a8 



{dki-KfdfdjT^djT^M^M^ 



3a6 

{dM7r)'dfdj7:dj7:MlMi 



3a4 6a8 (2M| + H^eMl> - SHMf) 6a^ (2M| + H^eMl> - 3HM^) 



9a6 



18a8 



36a8 



(75) 



where the above terms besides —{L^ + L4) are all at fourth order in perturbations. 



8. Symmetries 

Having written the complete Hamiltonian, we now proceed to calculate the four-point 
function contributions arising from interaction terms at third and fourth order. We 
employ here the IN-IN formalism [26, 27, 28, 29] and conveniently split the contributions 
to the four-point function as the ones arising from terms that make up the contact 
interaction diagram and the ones that generate the scalar exchange diagram as in the 
figure below. 
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Figure A: On the left, the scalar exchange diagram. Contact interaction diagram on the right 
It is useful at this stage to offer some comments on the calculations we are going 
to present. As mentioned, the literature already contains a thorough analysis of the 
trispectra for general single-field inflation models, see for example [53]. Work on the 
four-point function for ghost inflationary models has recently been presented [37, 38]. 
Our starting point, being based on a comprehensive effective theory, clearly encompasses 
all these models. Working with the effective Hamiltonian above translates into many 
immediate advantages as listed before but, on the other hand, in calculating the resulting 
four point function, one faces a substantial number of terms and it is therefore natural 
to look for some ordering principle which would single out some contributions to the 
trispectrum as the leading ones and allow us to concentrate on them only. In this context 
employing a symmetry for the whole theory can prove very useful. Indeed in [59, 38] 
the authors consider only those allowed by a particular (approximate in [59]) symmetry 
of the action, respectively: 

SI : TT -TT ; S2 : TT -TT and t ^ -t . (76) 

We plan here to employ our general effective theory to show that, allowing some freedom 
on the coefficients that modulate the various terms in the third and fourth order 
action, within each one of the two distinct and quite restrictive symmetry requirements 
above there are novel curvature-generated terms in the action that should not be disre- 
garded as negligible and that, furthermore, show some distinctive features in the shapes 
of the trispectrum. We will also describe terms allowed by both the symmetries in 
Eq. (76) combined. Of course, one need not employ symmetries to switch on or off any 
specific operator in the action. Most of the contributions are indeed freely adjustable 
by the correspondent M„ coefficient, a procedure which is, in principle, legitimate since 
the underlying theory is unknown. We choose here to restrict ourselves to considering 
only symmetry-abiding terms. Let us comment on each one of the symmetries. 
SI is built upon the following considerations. Often the same M„ coefficients multiply 
terms of different perturbative orders; consequently the amplitude /tvl of the 3-point 
function will be related to the amplitude of higher order correlators, notably to tjvl; the 
amplitude for the four-point function. Whenever the leading part of the trispectrum is 
generated by these types of M„'s one can estimate that for its effect to be observable 
ttvl has to be five orders of magnitude larger than /tvl [59], which leaves little room for 
feasible models. On the other hand, one quickly realizes those M„'s whose first term 
starts only at the fourth perturbative order (M4,Mio...in Eq. (74)) are not plagued by 
this problem. This then represents a natural way to obtain inflationary models which 
allow a large, detectable trispectrum untied to the interactions which make up the bis- 
pectrum (which might well be small now)*. Indeed, in [59] the authors investigate on 
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the size of all the interactions driven by the M4, Mn, Ni** coefficients in Eq. (74) and 
show that the leading interactions driven by these parameters are all consistent with 
the TT — )■ — TT prescription and are expected to give a comparable signal ***. By con- 
struction then, the terms in the interaction Hamiltonian that are going to contribute 
to the trispectrum and be consistent with the reasoning that inspired the SI symmetry 
are only some of the ones that will make up the contact interaction diagram, namely 
those whose lowest order interaction is already at fourth order. This limits us to the 
contributions regulated by the following coefficients: M4, M10..M15, N1..N5. 
S2 symmetry, on the other hand, does not prohibit third order interactions, indeed in 
[38] the interaction 7r(V7r)^ is considered and, by inspection of Eq. (13), one can see that 
also other terms are allowed, the one regulated by M5 and, notably, the Mq n^dijuY /a'^ 
term. The Mg-driven term is particularly interesting because its contribution to the 
bispectrum calculations of [22] generates an interesting flat shape. The scalar exchange 
diagram will then be built out of the third order S2-obeying terms in the action. In 
particular, inspired by previous findings, we are going to give a detailed account of the 
Mq contribution. 

If both SI and S2 are to be enforced one must also exclude from the hst of Sl-abiding 
interactions the ones multiphed by Mio, M13, M14, M15. A more clear picture of the 
situation concerning the various symmetries is presented in Table 1 below. 



Table 1 
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S2 
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X 
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X 


/ 


/ 


/ 


/ 


/ 





The Coefficients marked with " /" in correspondence of a given symmetry S are 
S-invariant, those marked with "X" violate the S symmetry. 



*One needs also to check that the interactions driven by coefficients that multiply also third 
order fluctuations do not become important in the form radiative corrections to the bispectrum. 
This check is done in [59] and ensures that loop corrections of those terms are not relevant. 
**In the same spirit of the analysis done in [22] for all curvature- generated terms at third 
order, the authors of [59] consider in the v2 of their paper some extrinsic-curvature terms 
generated at fourth order. They also comment on their importance in near de Sitter limit and 
their conclusions apply to our Mii,Ni parameters. 

***It would be interesting to understand to what kind of models, in terms of the fundamental 
scalar field, the simple resulting effective Lagrangian corresponds in this case. 
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Note that each M„ coefficient might multiply many interactions at each perturbative 
orders and therefore we mark the coefficient as invariant under a symmetry when all 
the leading interactions it multiplies are invariant under Si or S2. Determing the 
properties of the coefficients in the second row requires no effort, as one can easily verify 
these M^'s first appear in the action as multipliers of fourth-order terms. Things are 
less linear with the coefficients in the first row (except for M4) as they appear at fourth 
order both multiplying bare interaction terms and multiplying other coefficients as well 
as interaction terms (for an example of the latter case see the terms written explicitly 
in Eq. (75)). They also appear at third and some also at second order in perturbations. 
One then must carefully check that, given a particular coefficient M„, in none of the 
interactions it multiplies at any order the leading terms violate the symmetry. For 
M2, Ms, Mr,, Mq in the first row one can verify after some checks that these terms all 
parametrize indeed approximately invariant interactions upon requiring the coefficient 
M| to be much smaller than the typical M„ such as M2..MQ. This is because in the 
fourth-order Hamiltonian in Eq. (75) there are terms of the form 



which one then assumes to be subleading. We stress this point because it emerges clearly 
and naturally in the effective theory approach. 

9. Trispectrum (amplitudes and shapes) 

9.1. IN-IN Formalism 

We are going to employ the IN-IN formalism to calculate the foTir point function of 
curvature perturbation. The most general and compact expression for such a quantity 



(78) 



where T and T indicate respectively anti-time order and time order operations, |0) and 
1^2) stand for the vacuum of the free and interacting theory. 

Expanding both the exponentials in Eq. (78), we single out the first non vanishing terms 
that will contribute to the scalar exchange and contact interaction diagrams. 



oc 




X {M|, M^, M^, M|} X (S2 - violating interaction) 



(77) 



is: 



(^^iaia2Cfc3a4(t)i^^) = 



{0\f{t / d'xdtn3{x)}CklCk2Ck3Ck4{t)T{-l / d'x'dt"n3{x')}\o) 




+ {0\f{'- I I d^xdtd'x' dt"'H3{x)'H3{x)}CklCk2Ck3Ck4{t)\0) 
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+{0\CkiCk2Ck3Ck4{t)T{^-^ r r d'xdt d'x dtH^{x)H^{x)m 

^ J — oo J — oo 



to 



+ {0\T{t / d'xdtH4{x)}CkiCk2Ck3CkA{tm 

J — oo 

+{Q\CkiCk2Ck3CkA{t)T{-t I ° rfV dt"H4{x')}\0) + ... (79) 

J — oo 

where "Ha, Ti^ are the third and fourth-order Hamiltonian in the interaction picture. The 
latter two terms make up the contact interaction diagram, the rest is responsible for the 
scalar exchange. Let us also remind the reader that the gauge invariant observable 
( is, at first approximation, linearly related to the scalar vr via ( = —Hit. Also, 
already at this stage one can see that the result of the four point function is going to 
depend on six variables. All wavefunctions, once in Fourier space, depend only on the 
magnitude of their momenta. There are at most ten fields involved in the contractions, 
eight of which will always depend on the magnitude of the four external momenta 
{ki, k2, k^, k^j. We are left with one last contraction between two fields depending on 
the magnitude of one vector which, by construction, is going to be the sum of two 
external momenta. It turns out that, employing the overall momentum conservation, 
two variables are sufficient to describe any of these linear combinations, we choose 
(ki2 = 1^1 + ^2 1; ^14 = 1^1 + ^4!), giving a total of six variables. As clear from above, 
the Mg-driven third-order interaction we are going to consider further depends on scalar 
products between the various momenta but, as one can easily verify, these can all be 
fully specified by using the six variables introduced above. All the variables we will 
employ are represented in the figure below. 



k 



2 



Figure B: the regular tetrahedron described by the four external momenta and ki2,ki4 



In order to get a tetrahedron as the one in Fig B one must enforce the following 
inequalities: 

cos{a - /3)> cos(7) > cos{a + /3) (80) 

with 

^2 I l2 _ 1.2 h'^ -U h'^ - i-2 i-2 , iL2 _ p 

cosia) = ^^±l^^; cos{P)= 't'l ^ = " \l\ ■ (81) 

Z/C1/C14 /fC2fvi4 ZK1K2 
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Prom —1 < cos{a, /3,'j) < 1 one also obtains the usual triangles inequalities. Here we 
single out some of the inequalities which we are going to use in \A'hat follows: 

fc?2 + A;i4 < 4 ; ^Jl + kl-2k4 < ku < ^1 + A;| + 2k^. (82) 

In order to have a visual intuition and understanding of the result, once the calculation 
of the several contributions to the trispectrum is performed one needs to set up a number 
of configurations in which four out of the six variables are held fixed. Having more than 
one configuration also increases one's ability to distinguish the signatures of different 
interactions. Following [53], we adopt the set up described below: 

• Equilateral configuration: all the external momenta have the same magnitude 
k — ki — k2 — kz — k^] the two variables left are plotted as ki2/k,ku/k. 
Note that when plotting in this configuration we will use the first inequality in 
Eq. (82). Incidentally, this is the only configuration for which exact calculations for 
the trispectrum in ghost inflation have been presented (see [38]) so far. Note also 
that for the equilateral as well as for the other configurations, one conveniently plots 
the result of the calculations in Eq. (79) for any specific interaction term multiplied 
by a factor of Y[i=i ^f- done also because this factor is generally common to 
all the contributions and so removing it sharpens the differences between the plots 
of each interaction term. 

• Folded configuration: here one has ku — >■ as well as ki — k2 and ka — k^. The 
second and third inequahties in Eq. (82) must be enforced in this case. The variables 
ku and k^ are the ones plotted in this configuration. 

• Specialized planar limit configuration: in this case we have ki — k^ — ku as well as: 



k 



12 



kl+'^ (k2k, + ^Jm-kl){Akl-kl) y\'\ (83) 



2kl 

The variables plotted are going to be k2/ki and k^/ki. 

Near double squeezed limit configuration: the tetrahedron is now a planar 
quadrangle and k'^ = k^ = ki2- The region of interest is in particular the one 
for which kz, k^, ki2 where the following relation holds: 

. /h'^f — h'^ 4- i-2 _|_ i.2\ _ jL2 p _|_ jL2 jL2 _|_ p p _|_ jL2 h.2 _ p p _ ^.4 , jL2p 
■\y l^iy 1^12 ' ""3 ' ""4/ "'sl"'s2 ' "'12"'14 ' "'12"'4 ' "'14"'4 "'14"'3 ""4 ' ru-^rij^ 



ko^ 



V2k4 



(84) 



with 

A;^i = 2y(^^^7Tkr^kI)(^^^4^^k7"^^ k^s2 = Vihh + ks • k4)(A;3A;4 - kg • k4). 

(85) 

In this well the last two inequalities of Eq. (82) will be imposed on the 

variables ku/ki and k4/ki. Note that, only in this configuration, what one actually 
plots is the result of Eq. (79) times Y[i=i instead of Y[i=i ^f- This is once again 
done in order to better appreciate the difference among the many interaction terms. 
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We now consider the result for the scalar exchange contribution focusing in particular 
on an interaction term (the one proportional to Mq in Eq. 13 ) case which proved very 
interesting in plotting the shape of the bispectrum [22] as seen in Section 4- 
In all the calculations that follow we use a simplifying assumption which has been 
verified to hold for 3-point functions and is expected to hold for higher correlators as 
well [22]. Instead of using the generalized wavefunction which comprises the general 
single-field inflation solution, the Ghost inflation one, etc. as its simplified limits, we 
employ the usual Hankel function H^i2{kcsT) as a solution to the equation of motion 
for the quadratic action. The rationale for such a simplification is that, as one can 
readily verify, the main contribution to higher order correlators comes as usual from the 
horizon-crossing region and precisely in that region the behaviour of the general solution 
of Section 2 resembles very closely the one of the simpler specific DBI wavefunction. 
We elaborate further on this fact in Appendix C where some examples and comparisons 
of explicit calculations are provided. 

Before moving to the detailed analysis of the shape-functions for several interactions 
terms in various configurations, let us comment briefly on the amplitudes generically 
associated with these interaction terms. As noted before, building on the freedom on the 
M„ coefficients allowed by the theory and on the possibihty of employing a small speed 
of sound, <^ 1 (the same holds for the parameters which represent the generalization 
of Cg, i.e. ao and /3o), one can obtain large values for the amplitude associated to each one 
of the curvature-generated terms we are going to study. This has been quantitatively 
verified for all the terms of Eq. (13) in [22]. As an example, consider the Mg-driven 
fourthorder interaction term. To estimate the size of the amplitude associated with a 
given interaction term one considers its ratio with the quadratic terms of the theory at 
freezing [21]. Applying this prescription to our example one obtains: 



£2 M| 7/2^^2 m4 ~t 



(86) 



where the linear relation C, = —Htt has been used; taking ( ~ 10~^ gives a rough 
estimate of the size of the non linear corrections. A number of useful consistency checks 
for trispectrum calculations have been outlined in the literature (see first and sixth 
references in [19]). We leave the task of performing such checks to future work. 

9.2. Scalar exchange diagram 

Here we are going to consider the interaction term Mq 'k{dijT^Y /a^. Note that, as 
opposed to the Mg, Mg-regulated terms in Eq. (13) which give a fiat shape for the 
bispectrum much like the Mg-driven interaction, this term is actually invariant under 
the symmetry S2 while for it to be (approximately) invariant under SI one needs 
to require its Mg coefficient to be much smaller than M4 which would in turn make 
its signal undetectable. We now write more explicitly the contribution of the Mq- 
driven third order interaction to the scalar exchange diagram. For all the details of the 
calculation, including contractions, we refer the reader to Appendix B. Consider here 
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just one particular contraction of the fields, the sample contribution we are after looks 
like the following: 

all contractions 



IT 



rt^ 



(87) 



^* * * * 



^^1— ^fei2^'^i^fc2(^i)(ki ■ ks 



tl 



C?i2 — 7i"fci2 7rA:3Vrfc4(t2)(k3 ■ k4) 



(88) 



where an overall momentum conservation delta and a factor of (27r)'^ have been omitted 
for simplicity. 

In Fig. 23, 24 below we plot what one obtains by summing over all contractions, 
accounting for the symmetry factors of the vertices and plotting the result. 
For the sake of comparison we often make reference to the shapes obtained in [53] and 
[38]. In the former work the so called local trispectrum is also plotted and compared 
with findings for general single-field inflation models; in the latter one shapes for 
the trispectrum of ghost inflation are presented in the equilateral configuration only. 
As mentioned, our starting Lagrangian comprises both these inflationary models; we 
decided to concentrate on plotting the novel curvature-generated terms that are invariant 
under SI, S2 or both. 





Figure 23. The equilateral configuration sliape is presented on ttie left from the scalar 
exchange contribution of the Mg-driven interaction. The shape-function is different 
from the plots presented in [53]. On the other hand, once a necessary change of 
variables has been performed, it is qualitatively similar to the shape for the contact 
interaction diagram which arises from the ghost inflation (Vtt)^ interaction term in 
[38]. 

On the right we plotted our findings for the Mg-generated interaction in the folded 
configuration. It very much resembles the ones obtained in [53] for the scalar exchange 
diagrams from DBI-like terms, especially from the interaction 7r(V7r)^ . In all the 
pictures above and below ki has ben set equal to unity without loss of generality. 
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3.0 



0.0 



Figure 24. On the left the plot obtained for the A/g interaction in the planar 
configuration plot. Note here some of the interesting features: as ^2,^4 — > the 
shape function goes to zero, much like it happens for DBI-generated interactions. As 
^2,^4 — > 2 the shape function reaches values that are negative, albeit slightly so. 
Looking at the k2 = line we see that it is convex, rather than concave as found for 
other interaction types in [53] 

On the right the planar limit double squeezed configuration is plotted. In the region of 
interest, namely for fci2 — > 0, the shape function is non-zero, finite and negative; this 
again is different than what found in [53] for a variety of DBI-originated terms. 

Overall we see that performing the shape analysis for this S2-abiding term we are 
able to find some distinctive features with respect to the DBI-generated contributions. 
This is quite interesting also given the fact that the very same Mg-Hiodulated interac- 
tion term gives rise to a fiat shape for the three point function [22] which contributes 
to enlarge the allowed classification of bispectra shape-functions for single-field models 
of inflation (with Bunch-Davies vacuum). 



9.3. Contact interaction diagram 

We now turn to the calculation of various terms that contribute to the contact interaction 
diagrams. The terms driven by M3, M4 are DBI-generated and have been calculated in 
a number of papers, notably [19, 53]. M2 is found in both DBI and Ghost inflationary 
theories [53, 38]. If we are to preserve both symmetries then we need to focus on 
Mil, M12, N1..N5 as one can easily check from Table 1. 

Given any interaction term at fourth order, its contribution to the contact interaction 
diagram can be written as: 



(7rfc,7rfe7rfc37rfc,(t ^ O))^^^' = -2X„ {Q\T{i / d'x dt H,{x)txI^'k1^ti1^t:IS ^ 0)}|0) 

J —00 

(89) 

where Xm stands for taking the imaginary part; a delta enforcing momentum conserva- 
tion and unimportant numerical factors have been omitted. 
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We have applied the above formula to a number of fourth-order interaction terms provid- 
ing some examples of notable SI and S2-abiding terms, SI or S2 invariant contribu- 
tions and, finally, terms that do not respect any of the symmetries above. We start with 
the S2-invariant interaction term ~ M|(V7r)*^. This contribution is present in both DBI- 
like and Ghost inflationary models (e.g. in DBI one simply has = 1 — 2M|/Mpif). 
Although already written down in [38, 53], this term has not been plotted in all four 
configurations described in section 4 we employ here. This is because in DBI-like theo- 
ries it is expected to be subdomimant with respect to the M^n^ term. 

• 01 = 1/2 M| (9i7r)^ /a^ 

For the shape function of the operator (Vvr)^ we see that a number of interesting issues 
arise. First, the plot in the equilateral configuration does not resemble any of those 
plotted in [53] *. Then, as we mentioned in Fig. 24, in the double squeezed configuration 
the ki2 — )■ limit gives a non-zero finite shape function. This is important because, up 
to the results in [53], this limit was thought as very useful to distinguish the leading 
contributions coming from interactions at third order in perturbations from the ones at 
fourth order in fluctuations. 



2.0 




Figure 25. The equilateral configuration shape for the Oi operator is presented on 
the left. It is different from the results in [53] and very much resembles the plot we 
obtained for the scalar exchange calculation. Notice that upon performing a change 
of variables, it is basically identical to the (Vtt)^ interaction term plotted in [38] (We 
elaborate further on this point in Appendix C). 

On the right we plotted our findings for the Oi interaction term in the folded 
configuration. As it will be for the other interactions, this configurations provides 
no particularly distinctive features that would allow to single out the constributions 
from the different interaction operators. 

In fact, all of the terms contributing to the scalar exchange diagram in [53] give a shape 
function which in the ki2 — )■ is finite. On the contrary, the leading contact interaction 
diagram contributions analyzed in [53] do vanish in this limit. 
* See App. C for a detailed account of the plot of this term first done in [38] with different 

variables and results essentially identical to ours despite a simplyfing assumption on our part. 
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Note also that if one is to relax the assumption of a Bunch-Davies vacuum, the authors 
of [53] showed that this is not true anymore. 




Figure 26. The Oi interaction planar configuration shape on the left is not 
exceedingly different from the ones presented in [53]: it vanishes for k2,ki — ^ 0, it 
is peaked for k2 = 2 = k^. On the k2 = k^ line the shapefunction is convex, rather 
than concave as for the contact interaction term plotted in [53] . 

On the right we plotted the Oi interaction shape function in the planar limit double 
squeezed configuration. Here we immediately note an interesting feature: despite this 
being a contribution to the contact interaction diagram, in the k^ = ki2 — > it gives 
a finite, non zero shape function. We comment more on this fact in the text. 

Next, we continue keeping our attention focused on terms which are SI and S2 in- 
variant. These include, in terms of their free coefficient, Mn, M12, iVi-.iVs. Since they 
generate shape functions which are qualitatively very similar, we chose to plot just two 
representative terms in this list. 

• 02 = 1/6 Ml^ 772 (9f 7r)2 /a^ 

In the plots below we see that the Mn-driven term, generates shapes which are 
very similar to the ones plotted in [53] for the DBI-generated term ~ tt^. 
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Figure 27. The equilateral configuration shape is presented on the left for the O2 
operator. Because of the way the space derivatives are written in Fourier space there's 
no ki2 , ki4 dependence and so one gets a plateau. On the right our findings for O2 in 
the folded configuration. 




Figure 28. The O2 planar configuration shape on the left: it vanishes for fc2, ki 0, 
it is peaked for k2 = 2 = k^. On the ^2 = k/^ line the shapefunction is now concave, 
just like one would get for the tt* interaction. 

On the right we plotted the shape function in the planar limit double squeezed 
configuration associated to the O2 interaction term. 

We proceed with the other representative term: 
• O3 = 1/4! N^^l^:^ij^l^jk^i^k^'K /a^ 

The differences with respect to the Mn-driven interaction shapes are to be find 
in the first and the third configuration: in the first configuration they are due to the 
A;-dependence of the interaction, on the third configuration gives a plot similar to 
the one tuned by the M2 coefficient (see Fig. 26). 
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Figure 29. The equilateral configuration sliape for O3 is presented on tlie left. On 
tlie riglit our findings for tire O3 operator in the folded configuration. 





Figure 30. The planar configuration shape on the left for O3. On the right we plotted 
the shape function of the ©3 interaction term in the planar limit double squeezed 
configuration. 

We now turn our attention onto terms which violate one of the symmetries, SI in this 
case. Indeed, we analyze the interaction ~ M|M|7T^(9j7r)^ /a^ 



C4 = (M|M|)/(2M| + M|,e/72 _ 3MfH) x ^2 (a^Tr)^ /c 
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2.0.0 



Figure 31. The equilateral configuration shape for O4 is presented on tlie left. On 
ttie right our findings for the folded configuration associated to O4. Both are very 
similar to we obtained for the Mn in Fig. 27 and to what was found for the M4 -driven 
interaction in the literature [53] . 

For this interaction term we see the interesting feature presents itself in the fourth 
configuration where the ki2 — )■ hmit gives a finite shape function. 




2.0 0.0 



Figure 32. The planar configuration shape for O4 on the left. On the right we plotted 
the C'4-generated shape ftinction in the planar limit double squeezed configuration which 
gives again a non-zero and finite shape function even for fourth-order interactions such 
as the one under scrutiny here. 

We now proceed to plot our findings for one more term, precisely the leading fourth- 
order interaction term among the ones driven by Mg. It is clear that whenever this term 
gives a leading third order contribution (something one can achieve given the freedom 
on most M„'s ), it violates SI. S2 however, is preserved by the leading terms associated 
to this coefficient at third and fourth order as is clear from Eq. (75) and Table 1. 

.05 = l/6M6(5fc7r)2(9,,7r)Va' 
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1.0 



Figure 33. The equilateral configuration shape for the operator O5 is presented on 
the left: this shape for the plot function has not been seen before in the equilateral 
configuration. 

On the right our findings for the interaction term in the folded configuration. 



As one can see from Fig. 33, the plot in the equilateral configuration has no analogue 
in the shapes of [53, 38] for this configuration. It is somewhat reminiscent of the shape 
obtained for the vr^ of [53] but again, we stress it was obtained in a different configuration. 
The results plotted in Fig. 34 show once again that it is not safe in theories more general 
than DBI to attribuite to the planar limit double squeezed configuration the role to 
provide a distinctive signature in the ki2 — )■ limit that would enable one to distinguish 
between third and fourth-order interaction contributions (see the discussion in [53]). 




Figure 34. The planar configuration shape for O5 on the left. On the right we plotted 
the Os-generated shape ftinction in the planar limit double squeezed configuration 
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10. Conclusions 



In this thesis we have summarized the work done in [22, 23, 24]. The raison d'etre of 
these investigations was to study as general as possible a theory of single-field infiation 
and characterize the various infiationary mechanisms it encompasses mainly according 
to their non-Gaussian properties. 

To that aim, we employed the effective field theory approach of [21], which turned 
out to be a very powerful tool at our disposal. As detailed in Section 2, we first 
solved the equation of motion for the effective Lagrangian and obtained a wavefunction 
that interpolates between the classical solution of known inflationary models such as 
P{X, 0)-models and Ghost inflation. Our solution is actually more general, it does not 
just interpolate between known models. From the wavefunction one easily obtains the 
power spectrum of curvature perturbations, its tilt and its running. As expected, these 
quantities now depend on more variables than the usual three slow-roll parameters, 
there are actually five, generalized slow-roll parameters. We stress here that the two 
additional degrees of freedom are generated by extrinsic curvature terms in the quadratic 
Lagrangian; these can be significant in specific cases such as near de Sitter limit, or for 
small values of the generalized speed of sound or a combination thereof. 
The natural way to proceed is to tackle non-Gaussianities. We did so for the bispcctrum 
in Section 3 and 4- Here we showed that there exist a whole region of the parameters 
where extrinsic curvature-generated interaction terms play an important, possibly 
leading, role in determining the amplitude Jnl- Even more interestingly, a number 
of such operators independently generate a shape-function which peaks in the flat 
configuration: a feature which is quite uncommon for single-field models of infiation. 
Prompted by these findings, we also calculated the running of the contribution to /jvl 
given by one of such curvature operators, the one driven by the Mq coefficient: it turns 
out that the running depends, as it should, on linear combinations of the generalized 
slow-roll parameters but also on a parameter, eg, which is a function of the coefficient 
Me and its time derivative. We showed that, upon imposing suitable and mild bounds 
on the coefficients driving the various interactions, the running of f^L is dominated by 
Cq without spoiling any of the interesting features on the three and four-point function 
generated by Mg and without affecting the leading value of the power spectrum. In 
particular, we showed that the running ung can be such that ung ^ ^'(e, r), s). 
Mimicking the bispectrum analysis, and guided by the requirement of some additional 
symmetries on the action as an ordering principle [38, 59], in Section 7 and 8 we 
singled out the distinctive features of the trispectrum one obtains when considering 
extrinsic curvature-generated terms of a very general fourth-order Hamiltonian obtained 
in Section 6. It is important to note that all of these intereactions allow, by construction, 
for a large trispectrum. Some of them present features which also emerge in DBI- 
infiation and Ghost infiation [53, 38]. We have focused mainly on third and fourth-order 
interaction terms which have not been analyzed before and whose analysis reveals novel 
interesting effects. We were able to show that, unlike what happens in P{X, 0) models. 
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the analysis of the double squeezed configuration cannot give a clear cut clue as to what 
kind of signal comes from leading third-order terms as opposed to fourth-order terms 
in perturbations. We found many interactions generating a shape in the equilateral 
configuration that mimics the behaviour of the ghost interaction term ( i.e. (Vvr)^ ) 
shape function first plotted in [38], which is quite different from the shapes of the DBI 
model (we also extended the Ghost inflation plots of [38] to three other conflgurations) . 
Triggered by the fact several interactions generated an interesting flat shape for the 
bispectrum, we considered the effect of one of those terms (for consistency we chose 
again the one controlled by the Mq coefficient) for the trispectrum. We calculated and 
plotted the contributions of this term to the scalar exchange and contact interaction 
diagram: a shape-function which has not been found before emerged in the equilateral 
configuration for the contact interaction contribution. 

Both at cubic and quartic order in perturbations we again and again came accross the 
realization of the following feature: a shape function which in general single field infiation 
models is only obtained either by employing a linear combination of operators (as far as 
the Bispectrum is concerned) or relaxing the Bunch-Davies vacuum requirement for the 
theory, quite naturally (Bunch-Davies vacuum, no linear combinations) arises in more 
general setups as the one employed here. Furthermore, it does so when considering 
several and independent interaction terms. 

All the above results clearly point to at least two important facts: first, the effective field 
theory approach has proven very fruitful in addressing the need to describe infiationary 
models from a unifying general perspective. Second, in the quest for predictions 
on important cosmological observables it is crucial to consider the effect of extrinsic 
curvature-generated interaction terms in the inflationary action. 

Prom here one might proceed in several directions. It is true that extrinsic curvature 
terms have shown several interesting and distinctive features. Is there a UV-safe theory 
that comprises a phase of expansion of the universe described by higher derivative 
interactions? Something of this sort is described for example in recent literature [60, 61]. 
All the effective field theory machinery has been used here assuming a Bunch-Davies 
vacuum for the theory. One might well ask what would happen if we were to start from 
an excited state, much in the spirit of [30]. Another possibility is to expand on the 
effective approach employed here: removing the shift symmetry requirement, it should 
be possible to describe from a very general perspective the so called resonance models 
which have generated quite an interest* in the recent hterature [62] -[66]. 



* These studies have revealed a somewhat unexpected and intriguing fact: small periodic 
features in an inflationary model can have important consequences on the non-Gaussianities 
of such model. 
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11. Appendix A 



Explicit expression for the slow-roll parameter er 

We give here an explicit expression for the time dependence of the slow-roll parameters 
which we called er- This quantity if first written in terms of the M, M parameters: 



m = 7 



5 

4 + 



/ 4(j/2M2e+iMi4) 



X 



Mo4 



!(i?2M|,6+2M24+3Mi4) 



7 



5 

4 + 



where 7[t] is the Euler function. From here one calculates the quantity er = 
obtaining: 

/ / 



(90) 
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-1- 



M2e+4M2'*+6Mi4 



2H^Mj,e + 



X 



/ M ^ 

I (2H^Mle + Mi^) + 2 (i^^M^ e + 2M2^ + 3M/) . / — 

V P ^ 1 \ P ^ 2 ^ lyy 2mMle + 4M2^ + 6Mi 

V ^ dt V2//2M|.e + 4M24 + 6Mi4 



X 



+7i/2M|6Mi^- f , , + 2M2^Mi^- ( , , 

^ dt\2H^Mj>e + AM2^ + 6Miy dt \2H^Mj,e + AM2'^ + 6Mi 

( , , j] - IQHMleM^'H , , 

cii V '^H^Mle + 4M24 + GMi^ ^ 2H^Mle + 4M24 + GMi^ 

/\/f 4 71//- 4 

^ 2i/2M|,e + 4M24 + GMi^ 2i/2M2 e + 4M24 + GMi^ 
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M ^ 

i\/o'^ (^4 (H-J\If,f + 2ilf2^ + 3ilfi^) il'/o - Mo [hMj, {2cH + Hc^ + ^M^Mi + 123/i'^Jl/i)y 

2 {mmle + 2M24 + 3Mi4) 2 



[H'^Mle + 2M2^ + 3Mi^) ^ — + 16 [H'^Mle + 2M2^ + 3Mi^) ^ x 

^ ^ ^ ^ ^ 2Hm^e + 4M24 + 6Mi4 ^ p ^ 2 ^ 1 ; 



Mo^ 

2mMle + 4M24 + 6Mi4 

Explicit expression for the © functions of Eq. (35). 

IH^Mle + 8 (2M24 + SMi^) 
2 {H'^Mle + 2M24 + 3Mi4) ^ 
m^Ml^e^ - 16H^Mj,e (2M2^ + SMi^) + i/^jvfl, e {-UH^Mj,e + 8M2^ + ISMi^) _ 

2 {H^M^e + 2M2* + 3Mi4) (2H^M^e + Mi^ + Mq V2^^M|,e + 4M2^ + GMi^) ' 

4 {Hmi.e + 2M2^ + 3Mi4) ^ 

{m'^Mle + 24M2^ + SSMi^) i/^M^ e 

4 {H^Mj,e + 2M2^ + SM^^) (2H^Mj,e + M^^ + Mo'^y/2H^Ml,e + 4M2^ + GMi^) ' 

e 2M2^ 
^ e + 2M24 + 3Mi4) 

{m'^Mle + 3Mi4) 2M2^ 

{H'^Mle + 2M24 + 3Mi4) (2H^M^€ + Mi^ + MQ'^-s/2H^Mle + 4M24 + GMi^) ' 

^ e + 2M24 + 3Mi4) 

(+3Mi^ - AH'^Mle + 4M2^) 3Mi^ 
{H'^Mle + 2M24 + 3Mi4) (2H^Mj,€ + Mi^ + Mo^y/2H^M],€ + AM2^ + GMi^^ ' 

e.^, , (92) 

Explicit expression for variables in Eq. (38), (40). 
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AH (^-2M2* - 3Mi4 + M|,if^ 



(93) 



ao = ao(l + ^) ; ^0 = /3o(l + ^) ■ (94) 

In going from Eq. (38) to Eq. (40) it was safely assumed that a(T) ~ —jff^^z^, which 
in turn means that the cmplete and most general expression for the parameter xq is, to 
first order in generalized slow-roll parameters, given by: 

^ _ e ^ QM2^M2 ^ QMi^Mi 

3^^^ ^2e+'-e M 1.(95) 



AH(^-2M2^ + 3M^^ + Mj,H^ 2 2M| - ij + 3M4 
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12. Appendix B 

A more detailed presentation of the scalar exchange calculation for the Mg-driven term 
is presented here. We start from the scalar exchange part of Eq. (79): 

(^ICfeiCfc2Cfe3Cfc4(^)l^)s.e. — 

{0\f{i r Sxdt'H^{x)}CkMksCUt)T{-i r Sx'dtn^{x)m 

J —oo J —oo 



oo J — oo 



+(0|T{- y y df'xdtd'x dt Hs{x)H3{x)}Ck,Ck,CksCkAtm 

H^iCkMksCkMn^-^ f r d'xdtd^'x dtu3{x)U3{x)m. 

^ J —oo J —oo 

(96) 

Using Wick contraction on a generic operator 0, one has: 

T{0(ti)0(t2)0(t3)0(t4)} = N{(l){ti)(l){t2)(t){t3)(f){U) + all cofitractions} , (97) 

where N is the normal ordering operator. Note also that, being our vr ~ C operators 
squeezed between two vacua of the free theory, this reduces to considering only terms 
which are writeable as fully contracted contributions. For the anti-time order operator 
the same formula holds, only one needs to define contractions differently. We show 
below this difference: 



I 1 

0(fl,tl) 0(f2,t2)T = [<P^{Xl,h),(p'{x2,t2)]e{ti-t2) + [0+ (f 2 , ^2) , 0" (^1 , ^l)] ^(^2 " ^1 ) 

(98) 



I 1 

4){Xi,ti) 4>{x2,t2)f = [(p^{x2,t2),4>~{Xi,ti)]e{ti -ta) + [<P'^{Xl,ti),(j)~{x2,t2)]6{t2 - h) 

(99) 

where 

/d^k - r d^k 

j^Hk, t2) a^e^^-^^- r(x2, t2) = J j^,<l>*{k, i^) ale-''-'^ (100) 

Using the definitions above one gets several different contributions from Eq. (96). Note 
also that, using time and anti-time order definitions, the last two lines of Eq. (96) are 
just each other's conjugate and can therefore be grouped together. We now procede 
to write an explicit expression for the four point function generated by the Mg-driven 
contribution to the scalar exchange diagram: 

(7rfci7rfc27rfc37rfc4)l^g oc (101) 

-IT 4 ■ I 7r^i7r^2vrfe37rfe4(0) / (iti— 7rfci27rfci7rfc2(ki ■ k2)M dt2— 71112^3^4(^3 ' 
•J L V J-00 ^ J -00 ^ 

+ 5 permutations) + 
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I 2 ■ 2 ■ 7rfci7r^2^A:37rfc4(0) / (itl— 7rfel7rfel27rfe2(-ki2 • k2)^ / dt2—Tili2Trl37rl^{k3-k4 

+11 permutations) 

[ * * . 2 r ■* * * 2 

V J-00 « i-00 a 

+ ii permutations) 

4 • 7''fcl7rfc27rifc37rfe4(0) / dti^nklTTkl2T^k2{-k.i2 ■ 'is.2)'^ Cii2^7r^37rfei27r^4(ki2 • k4)^ 

J-00 a J-00 a 



permutations)] + 



2M^ 



3 

+ 5 permutations) + 



4 ■ ( 7rfci7rfc27rfe37rfe4(0) / rfii-j7rfei27rfei7rfe2(ki • k2)^ / C?i2-j7rfei27rfe3^fe4(k3 • k4)' 



I 2 • 2 • 7r^i7r^27rfe3^fc4(0) / C?^i-j7rfei7rfei27rfc2(-ki2 • k2)M dt2--jnki2'^k3'^k40^3 ' 
\ J-00 ^ J-00 ^ 

+11 permutations) 

2 • 2 • 7r^i7rfe2<37^fc4(0) / C?^i-j7rfei27rfci7rfe2(ki • k2)M (it2-j7i"fe37i"iki27rik4(ki2 • k4)^ 



+11 permutations) 

^3 /-ti ^3 



r « 2 / ' 2 

4 ■ 7rfei7rfe2Vrfe37rfc4(0) / dti-^TlklT^ll2T^k2{~ki2 ■ ^2)"^ \ (it2-j7rfc37rfcl27rM(kl2 • k4)^ 
J-00 ^ J-00 ^ 

+^5 permutations)] . (102) 

One then performs these calculations and plots the results to obtain Fig. 23,24. The 
situation for the contact interaction diagram contributions is considerably simpler as 
there is just one time intergral to be performed and two less fields to be taken into 
account. 
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13. Appendix C 



We want here to show with an example what seems to be a general feature concerning 
the use of (reasonably) approximated wavefunctions in the calculation of higher order 
correlators. In [22] we found that in performing an exact calculation for correlators 
in a very general theory such as the one we employed in this paper, whenever a 
given interaction term was producing a shape function for the trispectrum which one 
could qualitatively classify equilateral, so was the calculation performed with a 

simplified wavefunction. This is due to two independent reasons. First, we start from 
the realization that, precisely in the horizon-crossing region, which is where one expects 
the main contribution to any n-point to come from, the exact general wavefunction [22] 
and the usual one, H^^}^(cskr), which in these theories is an approximated solution, 
behave very similarly. Secondly, in [22] we concluded that most of the distinctive effects 
of the bispectrum where due not to the particular k-modes dependence of the result 
of the integrals like the one in Eq. (88), but on the fraction of that k-dependence that 
could be taken outside the integral, so on the part of the k-dependence not directly 
attached to the time behaviour of the wavefunction and which is common to the exact 
and approximated wavefunction. 

We now compare the trispectrum shapefunction of a ghost inflation interaction term, 
(Vtt)"', performed with the exact ghost solution in [38] with the results wc obtain 
employing the approximated DBI wavefunction, just what we used in obtaining all 
the shape functions presented here. 
The interaction reads: 



Its trispectrum shapefunction obtained through the approximated methods has been 
ploted in Fig. 25. In order to compare it with the exact calculation of [38] we need to 
change variables and turn to: 



C2 = h-k2; C3 = h-k; C^ = -l-C2-Cs; ^12 = ^2(1+^ = ^2 - C2 - C3) 

(104) 



We now show the plots obtained by performing this change of variable on our 
approximated result alongside the plot obtained with the exact ghost wavefunction 
taken directly from [38]. 




(103) 



i=i j=i 
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Figure 35. On the left the approximated result. The two shapefunctions are 
qualitatively identical barring an unimportant numerical coefficient due to a different 
normalization. 



14. Appendix D: IN-IN formalism 

See Section 9.1 
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